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ABSTRACT 


Rocket  attacks  on  civilian  and  military  targets,  from  both  Hezbollah  (South 
Lebanon)  and  Hammas  (Gaza  strip)  have  been  causing  a  major  operational  problem  for 
the  Israeli  Defense  Forces  for  over  two  decades.  In  recent  years,  U.S.  forces  are  facing 
similar  attacks  in  Afghanistan  and  Iraq  against  both  remote  military  outposts  and  in  the 
heart  of  Bagdad  (“Green  zone”).  The  insurgents  are  using  mortars  and  short  range 
rockets,  whose  launch  platforms  have  very  low  signature  prior  to  launch.  The  insurgents 
have  adopted  a  "shoot  and  scoot"  tactic  making  it  hard  to  detect  them  in  time  to  retaliate 
effectively.  In  this  thesis  we  present  a  new  analytic  probability  model,  that  addresses  this 
tactical  situation.  The  defender’s  decision  tradeoffs  are  explored  and  quantified.  A  new 
counter  mortar/rocket  tactic  is  suggested  and  explored  using  the  probability  model.  An 
extended  simulation  model  is  developed  to  explore  the  situation  when  the  defender  is 
using  a  sensor  that  is  subject  to  false  positive  detections. 
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EXECUTIVE  SUMMARY 


The  Operational  Problem 

Terrorists  and  insurgents  are  using  short  range  mortars  and  improvised  rockets 
against  U.S.  and  Israeli  Defense  Force  (IDF)  civilians  and  military  targets.  The  weapons’ 
low  signature  makes  them  hard  to  detect  before  the  actual  launch.  The  mortar  tube  or  the 
rocket  launcher  are  light  and  easily  transported,  and  may  be  disposable  after  the  firing. 
This  allows  the  attackers  to  use  a  “shoot  and  scoot”  tactic  making  it  hard  to  detect  them 
in  time  to  retaliate  effectively. 

The  Scenario 

We  consider  the  following  situation.  An  attacker  fires  a  salvo  of  several  rockets  at 
a  defender.  The  defender  uses  a  ground  sensor  to  detect  the  rocket  launches.  Following  a 
certain  number  of  rocket-launch  detections,  the  defender  launches  a  missile  aimed  at  the 
attacker.  The  attacker  stays  in  position  at  the  launching  location  until  he  has  fired  the  last 
rocket  in  the  salvo.  Soon  after  the  last  rocket  of  the  salvo  is  fired,  the  attacker  disappears 
from  the  launching  site  and  becomes  invulnerable  to  the  defender’s  missile.  Since  it  takes 
some  time  for  the  missile  to  reach  its  destination,  the  missile  may  reach  the  attacker’s 
launching  location  after  the  attacker  has  escaped,  in  which  case  the  defender  has  failed 
his  objective  to  kill  the  attacker.  Therefore  the  defender  has  a  tradeoff  between  launching 
his  missile  earlier  but  less  accurately  and,  launching  it  later  with  more  accuracy,  but  at 
the  risk  of  being  too  late. 

The  defender’s  firing  rule  determines  how  many  rocket  detections  need  to  be 
acquired  before  launching  the  missile. 
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The  Model 


The  analytic  model  is  described  by  the  following  equation: 


success 

N 

'sy'\ P {missile launched  on  kth  rocket }  •  P{missile  hit  on  time  \  k }  •  Pjattacker  is  killed  \  k  &  timely  hit } 

a-= 1 


Where  N  is  the  number  of  rockets  in  the  salvo.  In  the  analytic  model  we  assume  that  the 
ground  sensor  has  no  false  positive  detections.  Later  on  this  assumption  is  relaxed  in  a 
simulation  implemented  in  SimKit. 

Current  Tactic:  “Aim  first  -  Shoot  later” 

Typically,  to  counter  rocket,  artillery  or  mortar  fire  (henceforth  called  “rockets”), 
a  ground  sensor  (e.g.,  radar)  is  assigned  to  continuously  scan  a  “hot”  sector,  detect  enemy 
rockets,  and  locate  the  position  of  the  attacker.  To  allow  the  use  of  high  precision 
weapons  against  the  rocket  launchers,  and  in  order  to  reduce  collateral  damage,  the 
defender’s  ground  sensor  attempts  to  gather  as  much  information  as  possible  about  the 
location  of  the  launchers  before  the  defender  retaliates. 

The  ground  sensor  tries  to  track  the  trajectory  of  a  rocket  for  as  long  as  possible. 
The  estimated  trajectories  of  the  rockets  are  recorded  and  averaged  and,  based  on  these 
averages,  the  launcher  location  is  estimated.  The  decision  to  open  counter-fire  aimed  at 
the  launcher  may  depend  on  the  accuracy  of  the  location  estimation.  The  larger  the 
number  of  rockets  observed  by  the  ground  sensor,  the  larger  the  sample  size  of  the 
estimated  trajectories  and  thus  the  higher  the  accuracy  of  the  launcher  position.  However, 
while  this  tactic  of  “wait-and-observe”  increases  the  accuracy  of  the  location  estimate,  it 
consumes  valuable  time,  during  which  the  attacker  may  escape  and  thus  render  the 
counter-fire  useless.  A  firing  rule  in  this  tactic  will  call  for  a  missile  launch  after  as  many 
rocket  detections  as  possible. 

Suggested  Tactic:  “Shoot  first  -  Aim  later” 

To  mitigate  the  possibility  that  the  counter-attack  will  be  executed  too  late,  we 
propose  a  new  tactic.  This  tactic  may  require  a  new  type  of  weapon. 
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Consider  a  single  rocket  launcher.  Following  the  first  detected  rocket  in  the 
attacker’s  salvo,  the  defender’s  ground  sensor  obtains  an  initial  rough  estimate  of  the 
launcher’s  location  and  immediately  launches  a  precision  guided  missile  (PGM)  directed 
towards  that  estimated  location  of  the  attacker’s  launcher.  While  the  prompt  response  will 
save  precious  time,  this  is  not  enough;  accuracy  is  also  needed  for  effective  counter-fire. 
The  accuracy  will  be  achieved  by  updating  the  location  estimates  of  the  attacker,  while 
the  PGM  is  airborne,  based  on  detection  of  subsequent  rockets  in  the  salvo.  This  may  be 
achieved  either  by  sending  information  gathered  by  the  ground  sensor  to  the  missile  via 
uplink  or  using  a  second  sensor  onboard  the  missile.  A  firing  rule  in  this  tactics  will  call 
for  a  missile  launch  after  as  few  rocket  detections  as  possible.  A  firing  rule  is  defined  as 
the  number  of  detections  recorded  in  a  certain  time  period. 

Conclusions/Results 

A  model  for  evaluating  the  effect  of  various  firing  rules  has  been  developed  and 
the  conditions  for  the  suggested  tactic  success  are  presented.  The  suggested  tactic  has 
been  evaluated  and  explored  under  several  assumptions.  The  defenders  tradeoffs  between 
firing  rules  have  been  quantified  in  terms  of  the  probability  of  killing  the  attacker,  and  the 
rate  of  incidents  where  false  salvos  are  incorrectly  identified,  which  result  in  wasted 
missiles  and  possibly  collateral  damage. 
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I.  INTRODUCTION 


Citizens  in  Israel  have  been  under  the  threat  of  short-range  rockets  and  mortar  fire 
for  more  than  two  decades.  Hezbollah,  a  terror  organization  based  in  South  Lebanon,  has 
used  “Katyusha”  (107/122  mm)  rockets  to  target  towns  and  cities  on  Israel’s  northern 
border  since  the  80s.  These  attacks  have  been  an  ongoing  security  problem  of  varying 
intensity  over  time,  from  periods  of  sporadic  shootings,  to  periods  of  intense 
bombardments  lasting  days  and  weeks. 

This  form  of  terror  attacks  aimed  mostly  at  civilian  concentrations  reached  a  peak 
during  the  Second  Lebanon  war  (July  -  August  2006).  During  that  one-month  period, 
over  4000  rockets  were  fired  towards  Israel’s  northern  region  [1], 

Dealing  with  short  and  medium  range  indirect  fire  (such  as  Katyusha  rockets)  has 
always  been  a  great  challenge  for  the  Israeli  Defense  Force  (IDF).  This  type  of  weapon  is 
mounted  on  a  small  truck,  hauled  by  mules  or  even  carried  by  the  operator.  These 
weapons  have  very  low  signature  and  therefore  are  very  hard  to  detect  before  they  are 
launched.  Furthermore,  the  rocket  operators  have  adopted  a  “Shoot  and  scoot”  tactic; 
immediately  following  the  shooting  of  a  few  rockets  they  hurry  away  from  the  scene, 
often  leaving  behind  the  cheap,  disposable  launcher  (which  may  be  no  more  than  a  metal 
tube).  Therefore,  the  response  time  of  the  defender  must  be  very  short  to  effectively 
target  the  rocket  operator. 

In  recent  years  U.S.  forces  are  facing  similar  attacks  in  Afghanistan  and  Iraq 
where  insurgent  forces  fire  improvised  mortar  and  rockets  at  both  remote  military 
outposts  and  major  urban  areas  (e.g.  the  “Green  zone”  in  Bagdad)  causing  both  casualties 
and  disorder. 

In  this  thesis,  we  study  the  problem  of  counter-rocket  tactics  by  developing  a 
probability  model  that  captures  key  factors  affecting  this  combat  situation.  In  particular, 
we  propose  a  new  tactic  and  examine  it  using  the  probability  model. 
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The  model  developed  in  this  thesis  is  positioned  within  the  rich  research  area  of 
missile  defense.  In  their  book  [2]  Eckler  and  Burr  discuss  a  series  of  models  exploring 
offense  and  defense  strategies  in  allocating  ammunition  for  single  and  multiple  targets. 
Soland  [3]  explores  defense  tactics  that  try  to  minimize  the  fraction  of  a  target  destroyed 
by  the  attack  when  the  defender  has  shoot-look-shoot  capability  and  can  execute  several 
sequential  engagements  to  intercept  the  attack.  In  this  thesis  we  deal  with  a  special  case 
of  missile  defense  where  first,  the  defender  attacks  the  launching  platform  or  its  operator 
rather  than  the  incoming  missiles;  and  second,  the  defender  has  to  tradeoff  timeliness  (i.e. 
the  target  is  time  sensitive)  with  accuracy.  We  introduce  a  new  tactic  (which  may  call  for 
a  new  weapon  as  well)  that  tries  to  mitigate  the  defender’s  dilemma  “wait  or  shoot”  by 
using  information  about  the  attacker’s  location  which  becomes  available  after  the 
defender  has  counter  fired. 

Some  similarities  can  be  drawn  to  Ravid  [4]  where  two  alternatives  of  defense 
against  attacking  aircrafts  are  compared:  early  engagement  before  Bomb-Release-Line 
(BRL)  versus  engagement  with  higher  kill  probability  After  BRL.  Another  similar  work 
is  by  Sweat  [5]  who  models  a  duel  scenario  where  two  combatants  approach  each  other. 
The  fire  of  each  combatant  becomes  more  accurate  as  they  get  closer  to  each  other.  Both 
combatants  must  decide  when  to  shoot.  In  our  case  however  only  the  defender  has  such  a 
decision  to  make  where  as  the  attacker’s  strategy  is  predetermined  and  unaffected  by  the 
actions  of  the  defender. 

The  thesis  is  organized  as  follows.  In  Chapter  II  we  provide  the  reader  with 
background;  we  describe  the  current  counter  rocket  tactic  and  the  suggested  new  tactic; 
and  finally  we  describe  the  model  scenario.  In  Chapter  III,  we  develop  the  general 
analytic  model  and  its  variation  under  different  assumptions;  we  present  the  results  and 
discuss  the  defenders  dilemma  “wait  or  launch”  under  several  assumptions.  In  Chapter 
IV,  we  explore  how  false  positive  detections  affect  our  proposed  tactics.  For  this  purpose, 
we  develop  a  simulation.  We  present  the  results  and  perform  sensitivity  analysis  for 
some  of  our  assumptions.  Finally,  in  Chapter  V,  we  present  a  summary  and  conclusions. 
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II.  THE  COMBAT  SITUATION  AND  ENGAGEMENT  TACTICS 


In  this  chapter,  we  provide  background  and  introduce  the  reader  to  the  counter 
rocket/artillery  combat  situation  and  the  current  tactic  employed.  We  discuss  the 
disadvantages  of  the  current  tactic  and  suggest  a  new  tactic  to  mitigate  these 
disadvantages.  Finally,  we  describe  the  scenario  we  are  about  to  model  and  introduce  our 
Measures  of  Evaluation  (MOEs). 

A.  CURRENT  COMBAT  SITUATION  AND  TACTICS 

Typically,  to  counter  rocket,  artillery  or  mortar  fire  (henceforth  called  “rockets”), 
a  ground  based  sensor  (e.g.,  radar)  is  assigned  to  continuously  scan  a  “hot”  sector  of  the 
battlefield,  detect  enemy  rockets,  and  locate  the  positions  of  the  hostile  rocket  launchers. 
To  allow  the  use  of  high  precision  weapons  against  the  rocket  launchers,  and  in  order  to 
reduce  collateral  damage,  the  defender’s  ground  based  sensor  attempts  to  gather  as  much 
information  as  possible  about  the  location  of  the  launchers  before  the  defender  retaliates. 

Specifically,  following  detection,  the  ground  sensor  will  try  to  track  the  trajectory 
of  a  rocket  for  as  long  as  possible  until  it  must  stop  either  because,  it  is  assigned  to  track  a 
new  rocket  or,  because  the  trajectory  of  the  current  rocket  is  too  high  to  be  “seen”  by  the 
ground  sensor.  The  estimated  trajectories  of  the  rockets  are  recorded  and  averaged  and, 
based  on  these  averages,  the  launcher  location  is  estimated.  The  decision  to  open  counter¬ 
fire  aimed  at  the  launcher  may  depend  on  the  accuracy  of  the  location  estimation.  The 
larger  the  number  of  rockets  observed  by  the  ground  sensor,  the  larger  the  sample  size  of 
the  estimated  trajectories  and  thus  the  higher  the  accuracy  of  the  launcher  position. 
However,  while  this  tactic  of  “wait-and-observe”  increases  the  accuracy  of  the  location 
estimate,  it  consumes  valuable  time  during  which  the  operators  of  the  enemy’s  launchers 
may  escape  and  thus  render  the  counter-fire  useless.  This  tactical  situation,  which  is 
described  in  the  introduction,  represents  situations  such  as  countering  mobile  artillery  or 
short  range  rocket  fire  by  insurgency/terrorist  attack. 
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B.  SUGGESTED  TACTIC  AND  WEAPON 

To  mitigate  the  possibility  that  the  counter-attack  will  be  executed  too  late,  we 
propose  a  new  tactic,  which  is  described  next.  This  tactic  may  require  a  new  type  of 
weapon. 

Consider  a  single  rocket  launcher.  Following  the  first  detected  rocket  in  the 
attacker’s  salvo,  the  defender’s  ground  sensor  obtains  an  initial  rough  estimate  of  the 
launcher’s  location  and  immediately  launches  a  precision  guided  missile  (PGM)  directed 
towards  that  estimated  location  of  the  attacker’s  launcher.  While  the  prompt  response  will 
save  precious  time,  this  is  not  enough;  accuracy  is  also  needed  for  effective  counter-fire. 
The  accuracy  may  be  achieved  by  updating  the  location  estimates  of  the  attacker,  while 
the  PGM  is  airborne,  based  on  detection  of  subsequent  rockets  in  the  salvo.  There  are  two 
possible  methods  to  create  and  utilize  these  updates: 

1.  The  ground  sensor  observes  the  trajectories  of  subsequent  rockets  and 
generates  improved  updates  of  the  location  estimates  of  the  enemy’s 
launcher  that  are  transmitted  to  the  missile.  Based  on  these  updated 
estimates  the  trajectory  of  the  defender’s  PGM  may  be  updated  in-flight. 

2.  The  missile  is  equipped  with  an  EO/IR  sensor  which  can  detect  the 
launches  of  subsequent  rockets  in  the  salvo  while  airborne,  and  based  on 
these  cues  it  autonomously  updates  its  trajectory. 

Although  the  ground  sensor  may  have  more  energy  and  resources  than  a  missile 
sensor,  it  is  possible  that  Method  2  may  allow  for  higher  accuracy  than  achieved  by  the 
fixed  ground  sensor  in  Method  1  because,  as  the  missile  approaches  the  attacker’s 
location,  the  range  gets  shorter  and  thus  the  location  estimates  becomes  more  and  more 
accurate.  If  these  in-flight  updates,  either  Method  1  or  Method  2,  are  feasible,  they  may 
compensate  for  the  disadvantage  of  launching  the  missile  as  early  as  possible  based  only 
on  a  rough  estimate  of  the  attacker’s  location. 

In  this  thesis,  we  study  under  what  technical  and  operational  conditions  this  tactic 
is  superior  to  the  current  one  where  the  defender  waits  and  gathers  as  much  information 
as  possible  about  the  location  of  the  launchers  before  launching  the  PGM. 
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C.  SCENARIO 

We  consider  the  following  situation.  An  attacker  fires  a  salvo  of  several  rockets  at 
a  defender  from  a  mobile,  possibly  improvised,  rocket  launcher.  The  defender  uses  a 
ground  sensor  to  detect  the  rocket  launches.  We  assume  that  the  ground  sensor  may  only 
detect  the  rocket  at  the  moment  of  launch;  it  cannot  detect  a  rocket  once  airborne  and  it 
does  not  gather  any  more  information  from  the  rocket’s  trajectory. 

Once  a  rocket  launch  has  been  detected  by  the  ground  sensor,  the  defender 
immediately  launches  a  missile  aimed  at  the  attacker.  We  assume  that  the  defender 
launches  at  most  one  missile  per  rocket  salvo.  The  attacker  stays  in  position  at  the 
launching  location  until  he  has  fired  the  last  rocket  in  the  salvo.  Soon  after  the  last  rocket 
of  the  salvo  is  fired,  the  attacker  (launch  platform  or  just  the  operator  of  the  weapon) 
disappears  from  the  launching  site  and  becomes  invulnerable  to  the  defender’s  missile. 
Since  it  takes  some  time  for  the  missile  to  reach  its  destination,  it  is  possible  for  the 
missile  to  reach  the  attacker’s  launching  location  after  the  attacker  has  escaped,  in  which 
case  the  defender  has  failed  his  objective  to  kill  the  attacker. 

If  the  missile  is  launched  in  such  a  time  that  allows  the  missile  to  reach  the 
attacker’s  launching  location  while  the  attacker  is  there,  then  the  attacker  may  be  killed 
with  some  probability.  The  (conditional)  kill  probability  of  the  missile  may  or  may  not 
depend  on  how  early  it  was  launched  during  the  salvo.  If  the  missile  is  “dumb,”  meaning 
that  it  does  not  have  any  in-flight  guiding  mechanism  that  can  improve  its  accuracy  based 
on  new  information  obtained  during  its  flight,  then  this  probability  is  fixed;  more  rocket 
detections  will  have  no  effect  on  the  accuracy  of  the  missile.  If  the  missile  is  “smart”, 
meaning  that  it  does  have  in-flight  guiding  mechanism,  it  may  use  information  while 
airborne,  gathered  either  by  the  ground  sensor  (Method  1  above)  or  by  an  onboard  sensor 
(Method  2).  In  that  “smart”  case,  some  of  the  residual  rockets  in  the  salvo  may  be 
detected,  which  will  update  the  missile’s  trajectory  and  improve  its  effectiveness  (higher 
probability  of  kill). 
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In  the  absence  of  false  positive  detections,  there  is  no  doubt  that  the  suggested 
tactic  will  improve  the  probability  of  kill  for  a  “smart”  missile,  which  inevitably  must  be 
more  sophisticated  and  expensive  weapon.  However,  every  sensor  has  some  rate  of  false 
positive  detections.  Implementing  the  tactic  as  suggested  above,  where  the  defender  fires 
a  missile  once  a  rocket  has  been  detected  by  the  sensor,  will  most  likely  lead  to  waste  of 
missiles.  In  an  extreme  scenario,  where  the  attacker  is  not  firing  at  all  (a  cease  fire)  the 
rate  of  wasted  missiles  equals  the  rate  of  false  detections.  Besides  the  cost  of  wasted 
munitions,  firing  needlessly  will  increase  collateral  damage.  We  therefore  introduce  false 
detections  to  the  model  in  Chapter  4  and  accordingly  suggest  a  missile  launch  decision 
mechanism  to  reduce  the  rate  of  unwanted  missile  launches. 

Our  major  measures  of  effectiveness  in  this  thesis  are  therefore,  the  unconditional 
probability  of  killing  the  attacker,  Psucesss ;  and  the  probability  of  a  false  salvo  declaration 
(and  missile  launch)  or  rate  of  false  salvo  declarations. 
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III.  THE  ANALYTIC  MODEL  -  NO  FALSE  POSITIVES 

DETECTIONS 

In  this  chapter,  our  main  assumption  is  that  the  ground  sensor  has  no  false 
positive  detections.  This  simplifying  assumption  allows  us  to  use  analytic  tools  and 
explore  our  MOEs  through  their  explicit  expressions. 

Notation  Definitions  and  General  Assumptions 

We  make  the  following  assumptions  and  notation: 

•  The  ground  sensor  has  no  false  positive  detections  (we  shall  relax  this 
assumption  in  chapter  4). 

•  Number  of  rockets  A  in  a  salvo  is  fixed. 

•  The  number  of  residual  rockets  in  the  salvo  after  the  k,h  rocket  is  N-k. 

•  The  inter-firing  times  7^  between  rockets  in  a  salvo,  are  iid  with  cdf 
F(tf)  and  pdf  f(tf). 

•  The  escape  time  T,  is  the  time  that  passes  from  the  firing  of  the  last  rocket 
till  the  moment  the  attacker  disappears.  This  time  is  independent  of  Tf 
and  has  cdf  G(te )  and  pdf  g(te). 

•  The  missile  flight  time  rm  is  fixed. 

•  All  detections  occur  at  the  time  of  the  rockets’  launches. 

•  The  ground  sensor  detects  a  launch  with  probability  a  independently  from 
rocket  to  rocket. 

•  The  defender  launches  a  missile  immediately  following  the  first  detected 
rocket  launch. 

•  No  more  than  one  missile  is  launched  per  salvo. 
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•  A  timely  hit  occurs  when  the  missile  reaches  its  destination  before  the 
attacker  has  escaped. 

•  The  missile  probability  of  kill,  given  the  target  is  present  at  time  of  impact 
is  Pm  (and  zero  otherwise). 

A.  GROUND  SENSOR  AND  “DUMB”  MISSILE 

1.  Assumption 

In  this  subsection,  we  explore  the  model  where  the  defender’s  missile  is  “dumb”, 
meaning  that  once  the  missile  has  been  launched  it  can  neither  gather  nor  receive  from 
the  ground  sensor  any  more  information  about  its  target  location.  Therefore  it  cannot 
adjust  its  trajectory  while  airborne,  and  thus  cannot  enhance  its  accuracy.  In  our  model 
this  assumption  is  implicit  in  the  fact  that  Pkill  is  fixed  regardless  of  the  number  of  residual 
rockets. 

2.  Deriving  the  Model  Equation 

We  wish  to  compute  Psuccess  -  the  probability  that  the  attacker  has  been  killed. 

We  first  look  at  the  distribution  of  the  random  variable  TN  k  -  the  time  left  from 
the  detection  of  the  kth  rocket  till  the  attacker  has  escaped. 

Since  TN_k  is  the  sum  of  N-k+l  random  variables:  N-k  rocket  inter- firing  times 
and  the  escape  time,  the  probability  density  function  of  TN_k  denoted  hN_k(tN_k)  is  the 
convolution  of  (N-k+l)  density  functions,  N-k  density  functions  of  Tf ,  the  inter-firing 
times,  and  one  density  function  of  Te ,  the  escape  time. 

We  get: 

V,  (tN-k )  =  ((/*•••*/)*  g)(tN_k )  (0.1) 

N-k  times 

By  integrating  we  may  find  the  corresponding  cumulative  distribution  function 

H N-k  N-k  )  =  —  t N-k  } 
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The  probability  of  success  is  therefore  given  by: 

N 

Psucsses  =  X  C1  "  °  '  P(TN-k  *  T,n  }  '  PkUl  = 

k=\ 

N 

k  =  \ 


(0.2) 


3.  Deriving  HN_k(tN_k) 


To  facilitate  carrying  out  the  convolutions  in  the  following  derivation  I  will  use 
the  Laplace  transform  and  its  following  properties. 

•  Given  a  non  negative  function  / (x) ,  its  Laplace  transform  is  defined  as: 

oo 

<p(s)  =  L(f)  =  ^e~sxf(x)dx  (0.3) 

0 

•  The  Laplace  transform  of  a  convolution  of  two  functions  is  the 
multiplication  of  their  Laplace  transforms: 

L(f*g)  =  L(f).L(g)  (0.4) 

f*g  =  rl{L(f)-L(g)}  (0.5) 

•  The  Laplace  transfonn  of  the  integration  operator: 


(0-6) 

(0.7) 

(0.8) 


•  The  Laplace  transform  of  an  exponential  factor: 

L1{L(f)-(s-a)}  =  e-a‘f(t ) 

•  The  Laplace  transform  of  a  time  shift: 

L'{L(f)-ea,}  =  f(t-a)-U(t-a ) 

Where  U  ( t )  is  the  unit  step  function. 

3.1  First  Case  -  Uniform  Distributions 

In  this  case  we  assume  both  the  inter-firing  times  and  the  escape  time  have 
Uniform  distributions. 
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We  denote  the  mean  inter-firing  time  as  juf  and  the  mean  time  to  escape  as  //, . 
The  density  functions  and  their  respective  Laplace  transforms  are  then  given  by: 


=  i  2 a 


—  0  <  1^  -  /u  f  <  a 


<pf{s)  =  L(f )  =  —(e-^-a)s  -e{flf+a)s)  =  -e2as)) 

las  las 

—  0  <  It  -  u  \  <b 

g(te)  =  \lb  U  ^ 


cpg(s)  =  L(g)  =  =  J_e-^-b)s(l_e-2bS)) 

lbs  lbs 


(0.10) 


Substituting  (0.9)  and  (0.10)  into  (0.1)  and  using  the  relationship  in  (0.4)  we  get: 


L(hN_k)  =  (pf(sYN~k)  -(pg(s)  = 

r_L e~^f~a',s n  _ g-2as'> _  e-2fo) \ 

las  lbs 


(0.11) 


We  expand  the  expression  using  Newton’s  binomial  formula: 


L(hN-k)  = 


N~kJ  (: 2a)iN-k)s(N-k) 


N-k  N  —  k  1 

£(-!)"( - K2“~[e 

r?)  m  lbs 


iMe~b)s  „-(X+i)s 


_e-^e+D>S  j 


(0.12) 


Rearranging  we  get: 


£(V,)  = 


Z^1)" 


(lay  -lb]  (.9 


JN-k+l) 


N-k 


-[(jUf-a)(N-k)+2am+(jue-b)]s  ^-[( jUf-a)(N-k)+2am+(jUe+b)]s  - 


(0.13) 


We  notice  that  (0.13)  is  composed  of  three  parts  (separated  by  parentheses),  the 


first  part  is  some  constant  independent  of  s,  the  second  part  is- 


t(AT-jf-+l) 


and  the  third  part 


is  a  sum  of  exponents  dependent  on  s. 


It  can  be  shown  that:  L  ( 


w  1  , 


,(JV-i+l)  - 


(N-k)\ 
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Recalling  from  relationship  (0.8)  that  exponents  in  the  s  domain  are  translations  in 
the  time  domain  we  get  the  final  result: 

1 


N-k 


X(-d" 


{2afN~k)  •  lb 
rN-k ^ 


n= 0 


1 


(0.14) 


v  m  J 


(N  -k)\ 


in  \  (N-k ) 


Where: 


pm  =  ( juf  -  a)(N  -  k )  +  lam  +  ( jue  -  b ) 
P'l  =  (jUf  -  a)(N  -  k )  +  lam  +  (jue  +  b) 


HN.k(t)  = 


We  can  very  easily  carry  out  the  integration  and  get: 
1 


N-k 


Zev' 


n= 0 


(layN-k>  -lb 
N-k 


1 


(0.15) 


v  m  J 


(N-k  + 1)! 


• W-pm ) 


m  \(A^-A:+1) 


u(t-pm)-(t-p ;> 


m  \  (N—k+\) 


u(t-p :)] 


3.2  Second  Case  -  Normal  Distributions 

In  this  case  we  assume  both  the  inter-firing  times  and  the  escape  time  have  nonnal 
distributions.  We  denote  the  mean  inter-firing  time  as  juf  with  variance  aj  and  the  mean 

time  to  escape  as  jue  with  variance  cr .  In  this  case  TNk  has  a  normal  distribution  with 
mean  juN_k  =  (N  -  k)  ■  /jf  +  jue  and  variance  a2N_k  =(N- k)  ■  a]  +  a] . 

We  note  that  a  Normal  random  variable  may  assume  negative  values  which  are 
not  appropriate  in  the  case  of  time.  However,  it  can  be  easily  shown  that  for  the  scenario 
explored  in  this  paper,  the  negative  tail  of  the  normal  distribution  is  insignificant  and 
truncating  it  will  not  have  significant  effect  on  the  results. 

3.3  Third  Case  -  Exponential  Distributions 

In  this  case,  we  assume  both  the  inter-firing  times  and  the  escape  time  have 
exponential  distributions. 
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We  denote  the  mean  inter-firing  time  as  /  and  the  mean  time  to  escape  as  . 


;k 

(0.16) 

(0.17) 


,  h 

The  density  functions  and  their  respective  Laplace  transfonns  are  then  given  by: 

f(tf)  =  ’  ?/(J) = — 

Af  +  s 


g(0  =  ;  (pAs)=L(g)=  /L" 


\  +  s 

Substituting  (0.16)  and  (0.17)  into  (0.1)  and  using  the  relationship  in  (0.4)  we  get: 


L(hN_k)  =  <pf(sYN~k)  -(pg{s)  = 


X  (N~k) 


(Af+syNk)  Ae  +  s 


(0.18) 


We  define:  s'  =  Ae+s  and  substitute  into  1 1 ) 


L{hfi-k)  ~ 


A/N-k)-Ae 


A  (JV-*) 


K  1  (T/.  ~Ae) 


(N-k) 


{Af-Ae  +  s')(N-k)  s'  (Af-AjN-k)  s'  (Af-Ae+s')(N-k) 


A/N~k)-Ae 

(N-k) 


(A') 


<\(N-k) 


(0.19) 


(/i'+5') 


(N-k) 


We  notice  that  (0.19)  is  the  multiplication  of  three  expressions  (separated  by 
parentheses),  the  first  one  being  independent  of  s’,  the  second  one  is  just  the  inverse  of  s’ 
and  the  third  we  recognize  to  be  the  Laplace  transform  of  the  Erlang  density  function 
with  parameters  N-k,  and  A'  =  Af-Ae. 


We  recall  the  relationship  in  (0.7)  and  so  we  get  the  final  result. 


A/N~k)-Ae  N^-'  u  ,)t((Af-Ae)ty 

wo =*-*,/  w.yi- 1 ^(/ g)  y  .  ] 


(0.20) 


(Af-Aey-  — o  ml 

The  exponent  at  the  beginning  of  (0.20)  is  the  result  of  the  substitution  we  made 
in  (0.18)  and  the  relationship  in  (0.7). 


4.  Analysis 

Now  that  we  have  an  expression  for HN_k (tN_k)  we  can  calculate  the  probability 

for  a  timely  hit  and  therefore  -  by  factor  multiplication  -  the  probability  of  killing  the 
attacker. 
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We  would  like  to  explore  how  Psucsses  is  affected  by  variations  in  the  technological 

(e.g.  missile  speed,  ground  sensor’s  probability  of  rocket  detection,  missile  lethality)  and 
operational  (e.g.  defender  attacker  distance,  rockets  rate  of  fire)  parameters. 

Recall  the  expression  for  Psucsses . 

(3.2) 

k= 1 

Next  we  analyze  the  problem  with  respect  to  the  three  aforementioned  cases:, 
Uniform,  Nonnal  and  Exponential  distributions.  In  all  three  cases  we  make  the  same 
assumptions  as  described  in  the  introduction  to  Chapter  3  and  in  section  3.A.l.We 
analyze  the  effect  of  missile  flying  time  on  the  probability  of  success.  Recall  that  for  the 
“dumb”  missile  Pmcsses  is  proportional  to  Pm  therefore  we  can  focus  only  on  the 
probability  of  a  timely  hit. 

We  first  substitute  our  results  for  HN  k(tN  k  )  composed  of  uniform  distributions 
from  subsection  A.3.1. 

Uniform  Distributions 

The  scenario: 

There  are  N  =  5  rockets  in  the  salvo.  The  attacker  is  using  some  kind  of  short 
range  mortar  similar  to  the  U.S.  Army’s  M224.  The  inter-firing  time  is  uniformly 
distributed  between  4.5  to  5.5  seconds  (reasonable  for  this  kind  of  weapon  see  [6]).  The 
attacker  is  well  prepared  for  escape  prior  to  firing  the  salvo  and  therefore  the  time  to 
escape  is  uniformly  distributed  between  14  to  16  seconds;  a  =  0.7  ;Pm  =  1 
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Probability  of  success  for  different  missile  flight  times 


Figure  1.  Probability  of  success  for  different  missile  flight  times* 

*  Flight  times  are  in  seconds 

Figure  lpresents  the  probability  of  success  for  different  missile  flight  times.  We 
see  that  the  probability  of  success  is  a  monotone  non-increasing  function  of  the  missile 
flight  time.  This  comes  as  no  great  surprise  since  we  know  the  longer  it  takes  the  missile 
to  reach  the  target  the  higher  the  chances  are  that  the  target  will  escape  before  the  missile 
hits  the  ground. 

We  notice  that  the  decrease  is  occurring  in  5  steps.  The  last  drop  to  zero  occurs 
when  the  missile  flight  time  is  about  35  seconds. 

We  note  that  the  average  time  the  enemy  is  in  position  after  the  first  rocket  is  fired 
is  given  by:  (N  - 1)  •  juf  +  jue=  35  seconds. 
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To  understand  the  behavior  of  the  probability  of  success  as  shown  in  the  curve  of 
Figure  1,  let’s  explore  the  probability  for  a  timely  hit  when  there  are  N-k  residual  rockets, 
given  by  P{TN_k  >rJ  =  l-H N_k  (rm ) . 

The  following  Figure  presents  P{TN  k  >  rm)  for  different  numbers  of  residual 
rockets  N-k,  and  missile  flight  time  rm  .  The  scenario  parameters  are  as  in  Figure  1 


Probability  of  timely  hit  for  different  missile  flight  times  and  different  number  of  rockets  left 


Figure  2.  Probability  of  a  timely  hit  for  different  missile  flight  times  (in  seconds)  and 

number  of  residual  rockets* 

*  To  properly  read  this  figure,  print  in  color. 

We  notice  that  the  probability  of  a  timely  hit  given  the  number  of  rockets  left  in 
the  salvo  is  a  sharp  step  function,  in  particular  for  a  small  number  of  residual  rockets.  The 
drop  moves  in  the  positive  direction  of  the  rm  axis  for  increasing  number  of  residual 
rockets. 
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This  makes  sense  because  each  rocket  left  in  the  salvo  increases  the  average  time 
till  the  target  escapes  and  allows  more  time  for  the  missile  to  reach  its  target.  We  see  that 
each  one  of  the  curves  is  indeed  shifted  by  about  5  seconds  which  corresponds  to  the 
mean  inter-firing  time  of  the  rockets. 

We  go  back  now  to  analyze  the  behavior  of  the  probability  of  success  in  Figure 
IRecall  from  Equation  (0.2)  that  the  probability  of  success  is  expressed  by  a  sum  of 
arguments,  each  one  containing  the  probability  for  a  timely  hit  given  a  different  number 
of  residual  rockets.  Consider  now  the  region  in  Figure  1  for  values  greater  than  35 
seconds.  Looking  in  the  corresponding  region  in  Figure  2  we  see  that  the  probability  for  a 
timely  hit  is  zero  for  any  possible  number  of  residual  rockets  (0  to  N-l)  in  this  scenario. 
When  the  missile  flight  time  is  very  long,  even  detection  of  the  very  first  rocket  (leaving 
4  more  rockets  in  the  salvo)  does  not  leave  enough  time  for  the  missile  to  reach  the  target 
in  time. 

Now  consider  the  region  where  rm  is  between  about  30  and  35  seconds. 
According  to  our  scenario  the  average  time  between  the  first  rocket  launch  in  the  salvo 
and  the  moment  the  attacker  escapes  is  35  seconds.  On  average  there  are  30  seconds 
between  the  second  rocket  launch  in  the  salvo  and  the  moment  the  attacker  escapes,  and 
so  on.  Notice  in  Figure  2  that  for  this  [30,35]  region  of  rm  the  probability  for  a  timely  hit 

is  non  zero  if  there  are  4  rockets  left  in  the  salvo,  and  zero  or  close  to  zero  if  there  are  less 
than  4  rockets  left  in  the  salvo,  meaning  that  the  first  rocket  must  be  detected  in  order  to 
have  a  significant  probability  of  timely  hit.  In  Figure  1,  we  see  a  step  up  in  the  probability 
of  success  when  rm  <  35  seconds.  If  we  decrease  the  missile  flight  time  further  we  enter  a 

region  in  Figure  2  where  the  probability  of  timely  hit  is  significantly  greater  than  zero  if 
the  number  of  rockets  left  in  the  salvo  is  at  least  3.  This  means  that  in  the  expression  for 
the  probability  of  success  (0.2)  two  of  the  arguments  in  the  sum  are  now  non-zero.  This  is 
why  we  see  another  step  up  in  the  probability  of  success  as  we  decrease  rm. 

As  discussed  above,  the  probability  of  success  decreases  in  steps  when  increasing 
the  missile  flight  time.  However,  we  notice  that  the  size  of  each  step  is  not  the  same;  the 

last  decrement  (when  rm  increases  above  35  sec)  being  the  largest. 
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To  understand  why  this  happens  we  go  back  to  the  expression  of  the  probability 
of  success  in  Eq.(0.2).  We  notice  that  the  arguments  of  the  sum  are  multiplied  by  a 
discrete  function  of  k,  d(k)  =  (1  - a)k '  -  a-  Pkm .  Since  a  is  between  0  and  1  d(k )  is 
monotone  increasing  with  the  number  of  residual  rockets.  d(k)  controls  the  step  size  of 
the  curve  in  Figure  1  (probability  of  success),  we  can  now  see  why  the  returns  for 
decreasing  the  missile  flight  time  are  diminishing. 

The  initial  increase  in  the  probability  of  success  (when  the  missile  flight  time  is 
just  small  enough  to  give  some  chance  of  success),  is  given  by  d(k)  for  the  case  when 
k=l,  which  gives:  <7(1 )  =  a  ■  Pkill .  So,  if  reducing  the  missile  flight  time  is  difficult,  and  we 

want  to  maximize  the  initial  jump  in  the  probability  of  success,  we  should  increase  a  as 
much  as  possible. 

At  the  heart  of  this  qualitative  behavior  of  step  decreases  in  the  probability  of 
success  seen  in  Figure  1,  is  the  fact  that  we  chose  time  parameters  (inter-firing  and 
escape)  with  relatively  small  variances.  These  small  variances  lead  to  unifonn  cdf  close 
to  a  step  function.  The  uniform  cumulative  distribution  becomes  close  to  a  step  function 
whenever  we  reduce  the  range  between  its  two  parameters  relative  to  their  midpoint 
(mean).  This  occurs  because  the  pdf  will  become  close  to  an  impulse  function  and  thus 
the  cdf  will  become  close  to  a  step  function.  This  in  turn  makes  the 
distribution HN_k(tN_k)  resemble  a  step  function. 

When  we  increase  the  range  of  the  distribution,  we  increase  the  variance  relative 
to  the  mean  of  the  underlying  distributions  and  we  get  less  steep  CDFs.  This  will  change 
the  discrete  jump  nature  of  the  curve  we  see  in  Figure  1  and  make  the  probability  of 
success  plot  ’’smoother”. 

Figure  3presents  the  probability  of  a  timely  hit  for  the  uniform  distribution  case 
with  greater  variability.  Figure  4presents  the  resulting  probability  of  successes. 

The  scenario  is: 

There  are  N  =  5  rockets  in  the  salvo  as  before.  The  attacker  uses  the  same  weapon 
however  he  fires  at  a  less  regular  rate  and  the  inter-firing  time  is  uniformly  distributed 
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between  4  to  6  seconds.  The  time  to  escape  is  uniformly  distributed  between  10  to  20 
seconds.  Notice  that  the  means  of  the  underlying  distributions  are  still  5  and  15  seconds, 
respectively,  as  before,  a  =  0.7 ;Pm  =  1 


Figure  3. 


Probability  of  timely  hit  for  different  Missile  flight  times  and  different  number  of  rockets  left 


Missile  flight  time 


Probability  of  a  timely  hit  for  different  missile  flight  times  (in  seconds)  and 
number  of  residual  rockets  (increased  variability)  * 


*  To  properly  read  this  figure,  print  in  color. 
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Probability  of  success  for  different  missile  flight  times 


Figure  4.  Probability  of  success  for  different  missile  flight  times  (in  seconds) 

(increased  variability) 

Exploring  the  case  of  normal  distributions  gives  similar  results  to  the  uniform 
case.  Therefore  we  move  directly  to  the  exponential  case.  Results  for  the  normal  case  are 
presented  in  appendix  A. 

Exponential  Distributions 

An  interesting  case  is  the  case  where  both  the  escape  time  and  the  inter-firing 
times  are  exponential  random  variables.  In  this  case,  the  mean  and  the  variance  are 
coupled  and  the  ratio  between  the  mean  and  the  standard  deviation  is  always  1 .  Therefore 
we  never  achieve  the  qualitative  behavior  of  step  decreases  and  always  get  a  smooth 
curve  as  in  the  high  variability  case  above.  There  is  still  some  effect  of  diminishing 
returns,  however;  it  seems  the  “knee”  of  the  curve  is  only  achieved  for  a  relatively  low 
missile  flight  time.  Another  interesting  observation  is  that  for  the  exponential  case  the 
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probability  of  successes  is  quite  high  even  for  very  high  values  of  missile  flight  time. 
This  is  again  due  to  the  high  variability  imposed  by  the  exponential  distribution  because 
the  variance  is  the  square  of  the  mean. 

Figure  5and  Opresent  the  exponential  case.  The  mean  times  are  the  same  as  in 
Figure  3  and  Figure  4 


Probability  of  timely  hit  for  different  Missile  flight  times  and  different  number  of  rockets  left 


Figure  5.  Probability  of  a  timely  hit  or  different  missile  flight  times  (in  seconds)  and 

number  of  residual  rockets  (Exponential  case)  * 


*  To  properly  read  this  figure,  print  in  color. 
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Probability  of  success  for  different  missile  flight  times 


Figure  6.  Probability  of  success  for  different  missile  flight  times  (in  seconds) 

(Exponential  case) 

5.  Tradeoff:  Wait  or  Launch? 

Recall  that  the  probability  of  success  depends  on  two  factors:  timely  hit  and 
accuracy.  So  far  we  assumed  perfect  accuracy  and  only  looked  at  the  defender’s  attempt 
to  respond  as  quickly  as  possible  and  so  increase  the  probability  of  timely  hit.  Flowever, 
in  reality  accuracy  is  not  perfect  and  by  waiting  and  sampling  more  rockets  in  the  salvo 
the  defender  may  improve  his  accuracy  at  the  cost  of  reducing  the  probability  of  timely 
hit.  So  the  question  is  wait  and  improve  accuracy  while  reducing  the  probability  of  timely 
hit,  or  launch  with  lower  accuracy  but  higher  probability  of  timely  hit. 

Let  us  now  assume  that  the  defender  launches  the  missile  after  M  detections.  The 
missile  will  be  aimed  at  the  average  of  the  M  location  estimates.  The  serial  number  k  of 
the  rocket  after  which  the  missile  is  launched  has  a  negative  binomial  distribution  with 

21 


parameters  M  and  a.  We  further  assume  that  the  defender’s  missile  has  a  cookie  cutter 
damage  function  with  lethal  radius  R: 

[  p  r  <R 

d{r)  =  \Pn  (0-21) 

[0  o.co. 

and  that  the  defender’s  errors  of  the  estimates  of  attacker’s  location  after  each  rocket 

.2 

detection  are  iid  circular  normal  random  variables  with  mean  //  =  0  and  variance  <T  .  The 

_2 

_2 

n  —  o  .  A/  i  t 

average  error  of  M  estimates  will  then  have  a  mean  and  a  variance  M  .  We 

then  get  the  expression  for  the  probability  of  kill  (see  also  [7]): 

,  .  Mr2  _MR 2 

Pkiu  =  \\d(r)-p—^Q  r-dr-d<j)  =  p-(  1-e  2°2  )  (0.22) 

We  rewrite  Eq.  (0.2)  as  follows: 


y  _ _ (1-a) 


MR 2 

aM-{\-HN_k{Tm))-p\\-e2^) 


(0.23) 


0  presents  the  probability  of  success  for  different  missile  flight  times  and  different 
M  values.  The  scenario  is  the  same  as  in  Figure  1 ;  that  is,  uniform  inter- firing  and  escape 
distributions.  We  choose:  p  =  R  =  <j2  =  1 . 
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Figure  7.  Probability  of  success  for  different  missile  flight  times  (in  seconds)  and 

different  M  values.* 

*  To  properly  read  this  figure,  print  in  color. 

We  see  that  if  the  defender  chooses  to  shoot  as  early  as  possible  immediately  after 
the  first  detection  (M  =1)  then  even  a  slow  flying  missile,  which  takes  as  long  as  35 
seconds  to  reach  the  attacker,  may  still  hit  the  attacker  with  probability  of  -0.15.  For 
M  >  1  this  probability  is  0.  However  in  this  case  the  maximum  probability  of  success  is 
about  0.4  for  different  defender  missile  flight  times.  On  the  other  hand,  waiting  for  the 
second  detection  does  not  allow  the  defender  to  use  a  missile  which  takes  longer  than 
about  32  seconds  to  reach  the  target;  however  the  maximum  probability  of  success  is 
about  0.62. 
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We  see  that  there  is  a  tradeoff  between  the  defender  waiting  to  launch  until  more 
than  one  rocket  has  been  detected  and  thus  improving  the  probability  of  kill  vs.  the 
defender  launching  the  missile  as  early  as  possible  and  saving  valuable  time.  Waiting  too 
long  (M=5)  is  not  beneficial  because  the  loss  of  valuable  time  overrides  the  gain  in 
accuracy  and  the  net  result  is  a  low  probability  of  success.  The  dashed  black  line  marks 
the  maximal  value  over  all  5  curves.  We  see  that  if  the  defender  has  a  very  fast  missile 
(rm  <15  sec),  the  defender  should  wait  for  3  detections  before  launching  the  missile. 

However  if  the  defender  has  a  missile  which  takes  more  than  about  15  seconds  to  reach 
the  target  but  less  than  25  seconds,  he  should  wait  for  only  2  detections,  and  if  the 
missile  takes  longer  than  about  25  seconds  the  defender  should  launch  his  missile 
immediately  after  the  first  detection. 

There  is  another  element  influencing  the  probability  of  success  beside  the 
probability  of  timely  hit  and  accuracy.  This  other  element  is  the  probability  of  missile 
launch.  Notice  that  in  our  formulation  for  every  value  of  M  there  is  some  chance  that  no 
more  than  M- 1  rockets  will  be  detected,  in  which  case  the  missile  is  not  launched  at  all. 
Consider  the  region  in  Owhere  xm  =  0 .  Obviously  a  timely  hit  is  guaranteed  in  this  case. 
One  might  expect  the  curves  to  be  ordered  such  that  the  most  “accurate”  missile  ( M  =  5) 
will  be  on  top,  and  the  least  “accurate”  missile  (M  =  1)  at  the  bottom.  However  we  see 
that  the  curve  with  the  lowest  probability  of  success  is  the  one  where  M  -  5  and  the 
missile  is  the  most  accurate.  We  get  this  result  because  the  event  5  out  of  5  rockets 
detected  is  less  likely  than  for  instance  than  the  event  1  out  of  5  rockets  detected.  Figure  8 
presents  the  case  where  the  probability  of  rocket  detection  is  1 .  The  rest  of  the  scenario  is 
just  as  in  0  Notice  when  rm  =  0  the  curves  are  ordered  by  M  from  top  to  bottom.  We 

notice  that  in  this  case  the  defender’s  best  policy  changes,  and  for  low  values  of  rm  the 
defender  should  wait  for  4  or  even  5  detected  rockets  before  launching  his  missile. 
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Probability  of  success  for  diffrent  missle  flight  times  and  M  values 


20  25 

Missile  flight  time 


Figure  8.  Probability  of  success  for  different  missile  flight  times  (in  seconds)  and 

different  M  values.*  (Perfect  detection) 

*  To  properly  read  this  figure,  print  in  color. 

0  presents  the  probability  of  success  for  different  missile  flight  times  and  different 
M  values  for  the  exponential  case,  that  is  exponential  inter- firing  and  escape  distributions. 
The  mean  inter-firing  time  and  escape  time  are  5  and  15  seconds  respectively.  We 
choose:  p  =  R  =  cr2  =  1 
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Probability  of  success  for  diffrent  missle  flight  times  and  M  values 


Figure  9.  Probability  of  success  for  different  missile  flight  times  (in  seconds)  and 

different  M  values.*  (Exponential  case) 

*  To  properly  read  this  figure,  print  in  color. 

Just  as  in  Figure  6,  we  note  that  the  curves  in  the  exponential  case  are  smoother 
and  with  a  longer  “tail”  than  in  the  uniform  case.  This  is  due  to  the  greater  variability  in 
the  exponential  case  and  to  the  pdfs  shape.  Comparing  0  with  0  ,  we  find  great  similarity 
between  the  two  figures.  We  find  the  same  curves  intersect  with  each  other  in  both 
figures  (for  instance  M  - 1  intersect  M  =  2,3 &  4 in  both  figures),  a  lower  M  allows  the 
defender  to  use  a  slower  missile  and  so  on.  In  both  cases  we  see  that  given  a  very  fast 
missile  the  best  policy  of  the  defender  is  to  launch  after  M  —  3  detected  rockets;  for  a 
slower  missile  the  defender  should  launch  after  M  -  2  detected  rockets;  and  for  a  very 
slow  missile  the  defender  should  launch  at  M  =  1  .  However,  in  the  exponential  case  the 
region  of  rm  where  M  -  2  is  the  best  policy  extends  from  rm  ~  5  to  rm  ~  40  .  Even 
beyond  rm  ~  40  where  M  =  1  is  the  best  option  for  the  defender,  the  difference  in  the 
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probability  of  success  between  M  =  1  and  M  -  2  seems  insignificant.  It  would  seem  then 
that  in  the  exponential  case,  if  there  is  some  ambiguity  regarding  rm  the  defender’s  “best 
bet”  would  be  to  launch  after  M  =  2  detected  rockets. 

B.  GROUND  SENSOR  AND  “SMART”  MISSILE 
1.  Assumptions 

So  far,  we  have  assumed  that  once  the  missile  has  been  launched  it  cannot  get  any 
more  information  about  the  location  of  the  target.  In  our  model  this  assumption  was 
implicit  in  the  fact  that  PkiU  was  fixed  regardless  of  the  number  of  residual  rockets  or  the 

value  of  rm . 

We  now  assume  that  the  defender  has  a  “smart”  missile,  meaning  that  the  missile 
has  an  in-flight  guiding  mechanism.  The  “smart”  missile  may  use  new  information 
gathered  while  it  is  airborne.  This  new  information  may  be  gathered  either  by  the  ground 
sensor  or  by  an  onboard  sensor.  With  this  capability,  some  of  the  residual  rockets  in  the 
salvo  may  be  detected  and  be  used  to  improve  the  estimation  of  the  target  location.  The 
improved  estimate  may  enhance  the  accuracy  of  the  missile  and  result  in  higher  value 

of  Pm. 

Recall  that  for  the  “dumb”  missile  Pm  is  fixed  regardless  of  the  number  of 
rockets  that  are  launched  during  the  missile’s  flight.  The  question  now  is  what  will  be  the 
behavior  of  Pm  for  the  “smart”  missile. 

First,  we  notice  that  Pm  for  the  “smart”  missile  is  not  smaller  than  Pm  for  the 

“dumb”  missile.  This  is  true  because  the  “dumb”  missile  may  be  considered  as  a  special 
case  of  the  “smart”  missile.  In  this  special  case  no  new  information  is  gathered  (i.e.  the 
sensor  providing  the  new  in-flight  information  is  very  poor). 

Second,  we  note  that  Pm  must  be  monotone  non-decreasing  in  the  number  of 

residual  rockets.  The  larger  the  number  of  residual  rockets,  the  more  opportunities  there 
are  to  obtain  additional  observations  for  estimating  the  target  location  and  thus  improve 
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the  trajectory  of  the  defender’s  missile  and  enhance  its  accuracy  and  effectiveness.  We 
therefore  conclude  that  Pm  (N  -  k)  is  a  monotone  non-decreasing  function  of  N-k,  the 
number  of  residual  rockets. 

Finally,  we  note  that  PkjU  must  also  be  a  monotone  non-decreasing  function  ofrm  . 

Information  gathered  after  the  missile  has  reached  its  target  is  lost  and  cannot  change  the 
missile’s  trajectory.  Therefore  only  rockets  launched  during  the  missile’s  flight  can  help 
improve  its  effectiveness.  The  longer  rm  the  more  opportunities  there  are  for  detecting 
residual  rockets. 

To  capture  the  full  picture  we  should  consider  the  probability  of  kill  to  be  a 
function  of  both  the  number  of  residual  rockets  and  of  the  missile  flight  time.  However 
finding  an  analytically  tractable  closed  fonn  expression  for  Pkm(N  -k,rm)  is  difficult  (in 

most  cases);  we  therefore  adopt  a  different  approach,  using  a  simulation,  to  explore  our 
model.  The  simulation  is  described  in  detail  in  Chapter  IV. 


2,  Ground  Sensor  Guided  “Smart”  Missile  (Method  1) 


We  first  explore  the  case  where  the  “smart”  missile  receives  new  information 
gathered  by  the  ground  sensor.  Following  the  same  arguments  as  in  section  3.A.5  we  find 
an  expression  for  the  probability  of  kill  similar  to  the  one  in  Eq.  (0.22).  However,  now  we 
consider  not  only  the  number  of  rockets  detected  prior  to  the  defender’s  missile  launch 
but  also  the  number  L  of  rockets  detected  out  of  V  rockets  fired  while  the  missile  was 
airborne.  We  obtain  the  following  expression: 

(. M+L)R 2 

Pm=p\  1-e^M  (0.24) 


We  revise  (0.23)  and  obtain  the  following  expression: 


p 

success 


N  f  i  ,  \  N-k  V  f  \  (M+L)f? 

=  ZL_i  =  -{\-aTL-cf-p-{\-e  ^  ) (0.25) 

k=M\M  V  v=0  SoVV 


k=M\  J  v=0 

Finding  a  general  expression  for  p{v  =  v\  N-k,rm}  when  the  underlying 
distribution  are  Unifonn,  is  quite  difficult  and  so  we  will  use  the  simulation  as  mentioned 
above. 
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Let  us  look  first  at  the  case  of  the  smart  missile  launched  after  M  =  1  detected 
rocket.  Figure  lOpresents  the  estimated  probability  of  success  for  a  ground  sensor  guided 
“smart”  missile,  launched  after  M  =  1  detected  rocket.  Each  estimate  is  from  a  single 
simulation  run  containing  a  sequence  of  1000  rocket  salvos.  The  results  for  the  “dumb” 
missile  we  have  presented  in  section  3. A. 5  are  also  presented  here.  The  scenario  is  the 
same  as  in  Figure  1  (Uniform  interfering  and  escape  distributions).  We  assume 
p  =  R  =  <j2  =  1. 
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Probability  of  success  for  a  ground  sensor  guided  "smart"  missile  and, 
for  the  "dumb"  missile  launched  at  different  number  of  detections  (M) 


Figure  10.  Simulation  estimates*  for  the  Probability  of  success  for  a  ground  sensor 
guided  “smart”  missile  (M  =  1)  and,  calculated  Probability  of  success  for  the 
“dumb”  missile  launched  at  different  number  of  detections  ( M)  ** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color. 

We  see  that  the  probability  of  success  is  no  longer  a  monotone  decreasing 
function  in  the  missile  flight  time.  Instead,  it  has  a  mode,  starting  at  a  relatively  low  level, 
increasing  to  a  maximal  probability  of  about  0.8  and  then  dropping  to  zero. 

We  recall  that  in  this  scenario,  the  missile  has  to  be  in  flight  to  receive  any  new 
information  after  it  has  been  launched,  this  accounts  for  the  initial  increase  in  the 
probability  of  success  we  see  in  the  figure  above.  As  rm  increases  there  is  a  higher 
chance  that  additional  rockets  will  be  fired  during  the  missile’s  flight;  additional 
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observations  increase  Pm  and  finally  the  probability  of  success.  As  the  missile  flight  time 
increases,  eventually  the  probability  of  a  timely  hit  decreases;  this  in  turn  decreases  the 
probability  of  success  and  eventually  reduces  it  to  zero. 

Notice  that  the  “smart”  missile  with  M  =  1  is  always  superior  to  the  “dumb” 
missile  withM  =  1 .  This  is  because  it  will  always  have  at  least  as  much  information  as 
the  “dumb”  missile  when  launched  immediately  after  the  first  detection.  The  probability 
of  success  for  the  “smart”  missile  when  M  =  1  drops  to  zero  alongside  the  probability  of 
success  for  the  “dumb”  missile  when  M  =  1 .  When  M  =  1  the  “smart”  missile  is 
launched  as  early  as  possible;  its  chances  for  a  timely  hit  are  just  as  good  as  those  of  the 
“dumb”  missile  fired  immediately  after  the  first  detection.  For  values  of  rm  lower  than 
about  5,  the  probability  of  success  using  the  “smart”  missile  with  M  =  1  is  the  same  as 
that  for  the  “dumb”  missile  withM  =  1 .  This  makes  sense  because  when  rm  is  very  small, 

there  is  very  little  chance  for  another  rocket  to  be  detected  during  the  “smart”  missile’s 
flight.  The  “smart”  missile  will  then  be  aimed  only  based  on  the  single  rocket  detected  by 
the  ground  sensor  prior  to  the  missile  launch,  the  same  as  in  the  “dumb”  missile  case  with 
M  =1 . 

Recall  that  the  maximum  of  the  probability  of  success  for  the  “dumb”  missile 
over  all  values  of  M  ( max(Psa<xasj)  ),  provides  us  with  the  best  course  of  action  the 

M 

defender  can  take  using  the  “dumb”  missile.  Comparing  Psuccess  for  the  “smart”  missile 
with  M  =  1  to  max(f’i(ra,ss)  for  the  “dumb”  missile,  we  see  that  for  high  values  of  r  the 

M 

“smart”  missile  with  M  =  1  is  more  effective  than  the  “dumb”  missile,  regardless  the 
value  of  M.  However  below  a  certain  value  (about  10  sec)  of  rm  it  becomes  better  to  use  a 

“dumb”  missile  launched  after  M  -  3  detected  rockets,  than  using  the  “smart”  missile 
with  M  =  1 .  In  this  case  the  missile  is  effectively  so  fast  that  there  is  no  need  to  rush  and 
launch  it  early.  It  is  almost  guaranteed  that  a  very  fast  missile  will  reach  its  target  in  time. 
However  it  is  less  likely  that  during  its  short  flight  more  rockets  will  be  launched  and 
detected. 
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Notice  that  the  peak  of  the  probability  of  success  for  the  “smart”  missile  is  higher 
than  max(P  )  for  the  “dumb”  missile.  This  is  because  when  the  defender  uses  the 

M 

“dumb”  missile,  his  best  course  of  action  utilizes  the  information  from  at  most  3  rockets 
(he  may  choose  to  use  the  information  of  more  rockets,  but  the  advantage  of  more 
information  is  lost  to  the  disadvantage  of  launching  too  late);  however  the  “smart” 
missile  is  launched  as  early  as  possible,  and  if  it  is  in  the  air  long  enough  it  may  use  the 
information  of  possibly  up  to  N  rockets. 

We  see  that  firing  the  “smart”  missile  after  M  - 1  detected  rocket  is  not  always 
better  than  the  “dumb”  missile.  It  is  interesting  to  see  what  happens  if  the  defender  uses 
the  “smart”  ground  sensor  guided  missile  and  waits  to  launch  it  after  more  than 
M  - 1  detected  rockets.  Each  estimate  is  from  a  single  simulation  run  containing  a 
sequence  of  1000  rocket  salvos.  The  following  figure  presents  the  simulation  results  for 
the  probability  of  success  for  a  ground  sensor  guided  “smart”  missile,  for  different  values 
of  M.  Also  presented  is  the  curve  of  the  best  course  of  action  the  defender  can  take  using 
the  “dumb”  missile  we  have  seen  in  Figure  7  The  scenario  is  the  same  as  in  Figure  1 
(FIniform  inter-firing  and  escape  distributions).  We  assume  p  =  R  =  cr2  =  1. 
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Probability  of  success  for  diffrent  missle  flight  times  and  M  values 
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Figure  11.  Simulation  estimates*  for  the  Probability  of  success  for  a  ground  sensor 
guided  “smart”  missile  launched  at  different  number  of  detections  ( M)  ** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color. 

Comparing  Figure  7and  Figure  llWe  see  that  for  any  given  value  of  M  the 
“smart”  missile  is  superior  to  the  “dumb”  missile.  However  as  M  gets  higher  the 
difference  between  using  the  “smart”  vs.  the  “dumb”  missile  diminishes.  This  is  because 
the  “smart”  missile  improves  on  the  “dumb”  missile  by  taking  advantage  of  the 
opportunities  to  detect  residual  rockets  while  it’s  airborne.  As  M  increases  there  are  fewer 
such  opportunities.  In  fact  there  is  no  difference  between  the  “smart”  and  the  “dumb” 
missile  when  M  —  5  since  there  are  no  residual  rockets  for  the  “smart”  missile  to  draw 
on.  The  dash-dot  black  line  marks  the  maximal  value  over  all  5  curves  (for  the  “smart” 
missile).  This  line  provides  us  with  the  defender  best  policy  using  the  smart  missile.  We 
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see  just  as  with  the  dumb  missile  the  faster  the  missile  the  longer  the  defender  is  willing 
to  wait  (higher  M).  The  defender  should  never  wait  for  more  than  M  —  3  detected  rockets. 
There  is  no  advantage  in  using  the  smart  missile  when  xm  is  very  small  because  there  is 
little  chance  of  detecting  an  additional  rocket  launch  while  the  missile  is  airborne. 

We  notice  that  in  the  Uniform  distributions  case,  just  as  in  section  3. A. 5  that  the 
increase  and  decrease  in  the  probability  of  success  occur  in  steps  due  to  the  low 
variability  of  the  distributions  we  chose. 


When  the  underlying  distributions  are  exponential,  it  is  easier  to  find  an  explicit 
expression  for  equation  (0.25).  When  the  inter-firing  times  are  Exponential  iid  random 
variables  and  given  k  we  can  say  that,  V  -  min(  N  -  k,  X)  where  A  ~  Poisson(AfTm)  . 


Furthermore  L-min(N-k,Y)  where  Y  ~  Poisson(a  ■  Afrm)  is  a  Poisson  distribution. 


We  can  now  rewrite  (0.25)  and  get: 


N 


P 

success 


-I 

k=M 


k- 1 
M- 1 


N-k 


■  (1  -  a)k~M  ■  aM 


■(\-HN_k{Tm))-YJP{L  =  l\N-k,Tm}-{\-e 


{M+l)R 

2-2  )  (0.26) 


/=o 


Where:  P{L  =  l\N -k,rm}  = 


(Af-a-zJ 


N-k- 1 


1-  £  P{L  =  i\N-k,zJ 


0</<  N-k 

l  =  N-k 


Figure  12  presents  the  probability  of  success  for  a  ground  sensor  guided  “smart” 
missile,  alongside  the  results  for  the  “dumb”  missile  we  have  presented  in  section  3. A. 5. 
The  scenario  is  Exponential  inter-firing  and  escape  distributions  with  means  5  and  15 
seconds  respectively.  We  assume  p  =  R  =  a2  =  1 . 
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Probability  of  success  for  diffrent  missle  flight  times  and  M  values 
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Figure  12.  Probability  of  success  for  a  ground  sensor  guided  “smart”  missile  and,  for  the 
“dumb”  missile  launched  at  different  number  of  detections  (M)* 

(Exponential  case) 

*  To  properly  read  this  figure,  print  in  color. 

We  find  great  similarity  between  Figure  Hand  Figure  12The  ground  sensor 
guided  missile  is  always  superior  to  the  “dumb”  missile  launched  at  a  given  number  M  of 
detected  rockets.  For  higher  values  of  Mthe  “smart”  missile  loses  its  advantage  over  the 
“dumb”  missile.  However,  note  that  the  peak  probability  of  success  for  the  “smart” 
missile  has  shifted  from  to  lower  values  ofrm .  For  instance  for  M  =  1  ,  the  peak 

probability  of  success  has  moved  from  rm  ~  23  secs  (for  the  Uniform  case)  to  rm  ~  15  secs 

(for  the  Exponential  case).  The  peak  itself  is  lower  for  the  Exponential  case  (about  0.62) 
then  for  the  Uniform  case  (about  0.8).  Because  the  inter-firing  time  in  the  Exponential 
case  has  greater  variability,  it  may  assume  higher  values  than  in  the  uniform  case.  This 

35 


reduces  the  probability  of  detecting  residual  rocket  fired  while  the  missile  is  airborne. 
The  probability  of  timely  hit  is  lower  in  the  Exponential  case  than  in  the  Uniform  case 
(see  Figure  2and  Figure  5).  The  net  result  is  a  lower  peak  probability  of  success  than  in 
the  Unifonn  case. 


3.  Onboard  Sensor  “Smart”  Missile  (Method  2) 

We  now  explore  the  case  where  the  “smart”  missile  receives  new  information 
gathered  by  an  onboard  sensor. 

The  onboard  sensor  presents  us  with  a  new  situation.  We  now  have  to  take  into 
account  the  improvement  in  both  the  probability  of  rocket  detection  and  the  probability  of 
kill  as  the  sensor  moves  toward  the  attacker’s  location;  the  shorter  the  distance  between 
the  sensor  and  the  launching  site,  the  higher  the  detection  probability  and  the  accuracy  of 
the  target  location  estimate.  We  assume  that  the  onboard  sensor  is  able  to  detect  and 
measure  the  attacker’s  location  at  any  range  and  it  is  not  limited  by  its  footprint  or 
orientation. 


It  is  well  known  (see  [8], [9])  that  both  parameters  -  detection  probability  and 
accuracy  of  location  estimate  -  are  related  to  the  signal  to  noise  ratio  or  the  energy 
reaching  the  sensor  from  its  target.  We  denote  this  energy  as  E. 


We  take  the  general  relationship  between  the  probability  of  detection  and  E  from 
the  Albersheim  empiric  equation  (see  [8]).  We  will  therefore  use  the  following 
relationship: 


Pd(E)  = 


C,  •  exp(F) 


(0.27) 


1  +  Cj  •  exp(F) 

Where  Cj  is  some  normalization-factor,  and  Pd(E )  is  the  probability  of  rocket 


detection  when  the  sensor  receives  energy  E  from  the  launching  site. 


Fet  us  assume  that  the  errors  of  location  estimation  are  iid  circular  nonnal  random 
variables  with  mean  zero  and  variance  a2  (E) .  In  [9]  it  is  shown  that  the  variance  is 
related  to  E  as  follows: 
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(0.28) 


<x2(£)  =  — 

E 

Where  C2is  some  normalization-factor.  We  assume  that  the  sensor  will  use  a 
simple  average  to  integrate  the  information  from  several  rocket  detections. 

If  at  the  time  of  impact  the  sensor  has  detected  S  rockets,  each  with  energy 
Et  /  =  1  ...S  then  the  variance  of  the  averaged  location  estimation  will  be  given  by: 

a2  =  ^  cr2  (Ej ,) .  We  have  already  found  the  relationship  between  the  probability  of 

i 

kill  and  the  variance  of  the  averaged  measurements  in  (0.22).  Substituting  a2  and  (0.28) 
into  (0.22)  gives  us  the  new  expression: 


4// = /K1- exP(- 


R2 

2cr 


2  n  2 


2  n  2 


=  /?•(!- exp(- 


S-R 


2S§ 

i  hi 


-))  =  P-(l-exp(- 
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SZR 


2  2>!(£r) 


» 


(0.29) 


Finally  we  note  that  E  is  some  function  of  D  the  distance  between  the  sensor  and 
the  target.  The  relationship  E(D)U  may  be  appropriate  for  an  active  sensor  (e.g.  A 
radar  which  suffers  two-way  attenuation  losses).  A  passive  sensor  might  have  the 

D2 


relationship  E(D)U  — -  (e.g.  Electro-optical  sensor  which  suffers  one  way  attenuation 

T~\Z 


losses)  The  relationship  E(D )  may  represent  some  other  kind  of  sensor. 


Equipped  with  these  expressions  we  can  now  use  our  simulation  to  explore  the 
case  of  an  onboard  sensor.  Figure  13  and  Figure  13Figure  15present  the  single  detection 
probability  of  kill  (based  only  on  the  information  from  one  rocket  detection  from  that 
sensor)  and  the  probability  of  detection  for  different  types  of  sensors  as  a  function  of  their 
distance  from  the  target.  The  distance  is  normalized  such  that  D  - 1  is  the  distance 
between  the  ground  sensor  and  the  target.  We  assume  that  the  onboard  sensor  is  much 
smaller  than  the  ground  sensor  (due  to  limited  space  on  the  missile)  and  therefore  is 
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receiving  10  times  less  energy  than  the  ground  sensor  i.e.  Eonboardsemor(D  =  1)  =  s"""^ei,S01 

.  This  is  why  at  same  distance  from  the  target  the  ground  sensor  is  more  effective  than  the 
onboard  sensor.  The  normalization-factors  Cx  and  C2  in  equations  (0.27)  and  (0.29) 

respectively  are  chosen  such  that:  a\ Egroundsensor)  =  1  and  Pd( Egroundsensor)  =  0.7  ;  We 
assume  p  =  R  -  1 . 


Figure  13.  Single  detection  probability  of  kill  for  different  types  of  sensors 

Figure  14. 
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Figure  15.  Probability  of  rocket  detection  for  different  types  of  sensors 
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Figure  16  presents  the  probability  of  success  as  a  function  of  the  missile  flight 
time  for  different  types  of  onboard  sensors  whenM  =  1 .  The  case  of  the  ground  sensor 
guided  missile  whenM  =  1  is  presented  for  reference.  The  scenario  is  as  in  Figure 
l(Uniform  inter-firing  and  escape  distributions).  We  assume  p  =  R  =  cr2 {E ground sensor)  =  1 . 


Probability  of  sucess  for  diffrent  types  of  onboard  sensors 


Missile  flight  time 


Figure  16.  Simulation  estimates*  for  the  Probability  of  success  as  a  function  of  the 
missile  flight  time  (in  seconds)  for  different  types  of  sensors  (M  =  1)  .** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**  To  properly  read  this  figure,  print  in  color. 

We  notice  that  for  high  values  of  rm  all  types  of  onboard  sensors  give  lower 

probability  of  success  than  the  ground  sensor.  This  is  because  large  rm  implies  that  the 

missile  is  slow;  therefore  when  the  residual  rockets  are  launched  the  onboard  sensor  is 
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more  likely  to  be  relatively  far  away  from  its  destination.  Because  the  sensor  is  far  away 
from  the  target  the  overall  information  it  provides  (as  detennined  by  the  number  of 
residual  rockets  detected  and  the  measurement  errors  they  have)  is  less  than  the 
information  that  the  ground  sensor  would  have  provided.  However  as  rm  becomes 

smaller  (and  the  missile  is  faster),  we  see  that  the  probability  of  success  for  the  onboard 
sensor  increases  (although  not  in  a  monotone  fashion)  until  it  goes  above  the  probability 
of  success  for  the  ground  sensor  guided  missile.  Finally,  as  in  the  case  of  the  ground 
sensor  guided  missile,  when  xm  is  small  enough  the  missile  is  so  fast  that  no  rockets  are 

launched  or  detected  during  the  missile’s  flight;  the  resulting  probability  of  success  has 
the  same  value  as  if  the  missile  was  “dumb”. 

We  note  that  the  probability  of  success  now  has  more  than  one  mode.  For  instance 
for  all  onboard  sensors,  there  is  a  gradual  decrease  in  the  probability  of  success 
when  5  <  rm  <10,  then  around  rm=  10  there  is  a  sharp  increase  in  the  probability  of 

success.  As  we  increase  rm  from  5  to  10,  the  missile  becomes  slower;  the  rockets 
launched  during  its  flight  are  detected  further  away  from  the  missile.  For  this  reason  the 
rockets  become  harder  for  the  onboard  sensor  to  detect  and  when  detected  they  provide 
less  information  and  have  a  less  effect  on  the  probability  of  kill.  However  since  the  inter¬ 
firing  time  is  unifonn  between  4.5  and  5.5  seconds,  during  a  flight  time  of  less  than  10 
seconds  no  more  than  2  rockets  may  be  launched  and  detected.  As  we  increase  rm  beyond 

10  seconds,  the  opportunity  to  detect  a  third  rocket  arises;  the  extra  rocket  will  provide 
more  information  which  causes  a  spike  in  the  probability  of  success.  This  ripple  effect 
would  not  have  been  so  distinct  had  we  chosen  inter-firing  times  with  a  greater  range. 

For  the  scenario  described  above,  we  see  it  is  not  beneficial  to  use  an  on-board 
sensor,  unless  we  have  a  relatively  fast  missile. 

Figure  17  presents  the  probability  of  success  as  a  function  of  the  missile  flight 
time  for  different  types  of  onboard  sensors  whenM  =  1,  for  the  Exponential  case.  The 
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case  of  the  ground  sensor  guided  missile  whenM  =  1  is  presented  for  reference.  The 
scenario  is  Exponential  inter-firing  and  escape  distributions  with  means  5  and  1 5  seconds 
respectively.  We  assume  p  =  R  =  cr2(Egl.oundsensor)  =  1. 


Probability  of  success  for  diffrent  types  of  onboard  sensors 


Figure  17.  Simulation  estimates*  for  the  Probability  of  success  as  a  function  of  the 
missile  flight  time  (in  seconds)  for  different  types  of  sensors  (M  =  1)  .** 

(Exponential  case) 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color. 

Notice  that  in  the  exponential  case,  when  the  onboard  sensor  receives  10  times 
less  energy  than  the  ground  sensor,  there  is  hardly  any  difference  between  the  probability 
of  success  if  the  missile  is  guided  by  the  ground  sensor  or  by  the  “Active”  onboard 
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sensor.  The  active  sensor  does  do  slightly  better  than  the  ground  sensor  for  rm  ~  10  secs. 
It  seems  that  in  this  case  investing  in  the  development  of  an  onboard  sensor  is  not 
beneficial. 

Finally,  we  explore  the  defender’s  option  to  wait  for  M  detected  rockets  before 
launching  the  missile.  The  next  4  figures  present  the  probability  of  success  for  the 
onboard  sensor  “smart”  missile  after  different  number  M  of  detected  rockets,  for  “Active” 
and  “Passive”  sensor.  The  scenario  in  Oand  Ois  Uniform  distributions  as  in  Figure  1  The 
scenario  in  Figure  20and  Figure  21  is  Exponential  inter- firing  and  escape  time  with  means 
5  seconds  and  15  seconds  respectively. 


Probability  of  success  for  onboard  sensor  and  diffrent  values  of  M 
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Figure  18.  Simulation  estimates*  for  the  Probability  of  success  as  a  function  of  the 

missile  flight  time  (in  seconds)  for  onboard  “Active  sensor,  for  different  values  of 

M.**  (Uniform  case) 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000  rocket 
salvos. 

**  To  properly  read  this  figure,  print  in  color 
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Probability  of  success  for  onboard  sensors  and  diffrent  values  of  M 
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Figure  19.  Simulation  estimates* **  for  Probability  of  success  as  a  function  of  the  missile 
flight  time  (in  seconds)  for  onboard  “Passive  sensor,  for  different  values  of  M* 

(Uniform  case) 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000  rocket 
salvos. 

**To  properly  read  this  figure,  print  in  color 
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Probability  of  success  for  onboard  sensor  and  diffrent  values  of  M 
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Figure  20.  Simulation  estimates*  for  the  Probability  of  success  as  a  function  of  the 
missile  flight  time  for  onboard  “Active  sensor,  for  different  values  of  M.** 

(Exponential  case) 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000  rocket 
salvos. 

**To  properly  read  this  figure,  print  in  color 
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Figure  2 1 .  Simulation  estimates*  for  the  Probability  of  success  as  a  function  of  the 
missile  flight  time  for  onboard  “Passive  sensor,  for  different  values  of  M.** 

(Exponential  case) 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000  rocket 
salvos. 

**To  properly  read  this  figure,  print  in  color 


Just  as  with  the  “dumb”  and  the  ground  sensor  guided  missile,  the  defender 
should  never  wait  for  more  than  M  =  3  detected  rockets  in  this  scenario.  The  black 
dotted  curve  is  the  “best  policy”  curve  for  the  onboard  sensor.  Also  plotted  for  reference 
are  the  “best  policy”  curves  for  the  “dumb”  and  the  ground  sensor  guided  missiles.  Note 
that  the  passive  onboard  sensor  missile  is  always  inferior  to  the  ground  sensor  guided 
missile  (under  these  scenarios).  The  “Active”  onboard  sensor  missile  does  a  little  better 
than  the  ground  sensor  guided  missile  in  the  Uniform  case,  but  in  the  Exponential  case 
there  is  no  distinguishable  difference  between  the  two. 
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IV.  THE  CASE  OF  FALSE  POSITIVE  DETECTIONS  -  A 

SIMULATION  MODEL 

A.  INTRODUCING  FALSE  DETECTIONS 

1.  Assumptions 

In  Chapter  III,  we  assumed  that  the  ground  based  sensor  has  no  false  positive 
detections.  In  this  chapter  we  explore  how  false  positive  detections  (henceforth  called,  in 
short,  false  detections)  affect  our  proposed  tactics. 

Generally  speaking,  every  sensor  has  some  false  detections,  which  may  be  due  to 
processes  internal  to  the  sensor,  or  due  to  external  conditions  such  as  weather  or 
interference  caused  by  other  systems.  In  our  model  we  assume  that  the  times  between 
two  consecutive  false  detections  are  iid  random  variables  independent  of  whether  a  rocket 
has  been  fired. 

We  assume  that  only  the  ground  sensor  may  have  false  detections.  The  missile’s 
onboard  sensor  does  not  have  false  detections. 

We  assume  that  the  ground  sensor  makes  a  measurement  of  the  attacker’s  location 
with  circular  normal  distributed  errors  and  variance  a2 ,  even  when  the  detection  is  false. 
We  make  this  assumption  to  simplify  the  simulation.  Obviously,  in  a  real-world  situation 
false  detections  may  reduce  the  missile’s  accuracy.  However  this  assumption  still  makes 
some  sense.  In  a  real  situation  if  the  false  detection  is  extremely  far  away  from  the 
attacker’s  location  it  will  not  be  associated  to  the  salvo.  Furthermore,  the  simulation  does 
handle  the  cases  when  false  detections  trigger  a  missile  launch  while  there  is  no  ongoing 
salvo  at  all.  As  in  Chapter  III  ,  we  assume  that  the  defender  launches  at  most  one  missile 
per  salvo. 

2.  The  Firing  Rule  (M,f) 

In  Chapter  III,  we  assume  that  the  defender  fires  back  a  missile  as  soon  as  he 
detects  a  rocket  within  the  salvo.  In  the  presence  of  false  detections,  such  a  policy  may 
lead  to  unnecessary  missile  launches  whenever  a  false  detection  occurs.  Such  missile 
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launches  are  disadvantageous  for  two  reasons.  First,  they  cause  waste  of  missiles,  which 
may  be  limited  or  expensive;  and  second,  it  may  lead  to  unnecessary  collateral  damage  as 
missiles  may  hit  innocent  targets. 

To  mitigate  the  problem  of  unnecessary  missile  launches,  we  introduce  a  new 
response  policy  for  the  defender.  We  assume  that  the  defender  knows  that  the  attacker 
fires  his  rockets  in  salvos  of  several  rockets  at  a  time.  The  defender  takes  advantage  of 
this  fact  and  estimates  when  a  salvo  has  begun  by  looking  for  a  series  of  detections  close 
in  time.  Following  such  a  cluster  of  detections,  the  defender  makes  a  salvo  declaration 
and  fires  a  missile  in  response.  Such  a  policy  may  reduce  the  rate  of  launching  missiles 
unnecessarily  since  now  only  false  detections  which  occur  in  close  proximity  will  lead  to 
a  false  salvo  declaration  and  unnecessary  missile  launches. 

The  firing  rule  is  defined  by  the  pair  ( M  ,T )  as  follows: 

Looking  at  a  sliding  window  in  time  of  size  T  time  units,  the  defender 
declares  a  salvo  the  moment  the  sliding  window  contains  at  least  M  detections, 
counting  from  the  last  detection  accounted  for  by  the  previous  salvo  declaration. 

Detections  from  or  before  the  previous  salvo  declaration  are  not  counted  to  reduce 
the  number  of  salvo  declarations  during  the  same  salvo. 

We  define  a  positive  salvo  declaration  as  a  salvo  declaration  which  is  made 
during  a  real  salvo  (from  the  moment  the  attacker  is  in  position  until  he  has  fired  the  last 
rocket). 


3.  Problem  Definition 

For  M  large  enough  and  T  small  enough  the  firing  rule  may  reduce  the  number 
of  false  declarations,  however  it  introduces  two  new  problems.  First,  the  probability  of  a 
positive  salvo  declaration  is  decreased  and  second,  valuable  time  may  be  lost  waiting  for 
enough  detections  to  meet  the  firing  rule  criteria,  thus  decreasing  the  probability  of  timely 
detection. 
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We  therefore  want  to  find  the  firing  rule  pair  ( M,T )  which  gives  the  best  tradeoff 
between  low  false  salvo  declaration  rate  and  high  probability  of  success. 

B.  SIMULATION  DESCRIPTION 

Initial  attempts  to  describe  the  problem  in  the  presence  of  false  detections 
analytically  did  not  lead  to  tractable  closed  form  solutions  which  lend  insights  easily.  We 
therefore  chose  to  approach  the  problem  using  discrete  event  simulation. 

1.  Events  Graphs 

To  write  our  simulation  we  have  used  discrete  event  simulation  (DES) 
methodology  and  the  SimKit  programming  toolkit  (see  [10]).  DES  simulations  maintain  a 
list  of  scheduled  events.  Each  event  may  affect  system  variables,  schedule  future  events 
or,  cancel  future  events.  DES  simulations  are  represented  by  events  graphs.  The 
following  figure  presents  a  simple  example  of  a  DES  event  graph. 


Figure  22.  Event  graph  example 

In  the  simple  event  graph  above,  there  are  two  types  of  events  represented  by 
circles:  doFoo  and  doBar(a). 


The  description  of  each  event  is  described  in  a  rectangular  frame  attached  to  the 

event  by  a  broken  line.  doFoo  increments  the  system  variable  numberOfFoos  and 
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instantiates  the  variable  “a”.  doFoo  adds  “a”  (received  from  its  argument)  to  the  system 
variable  numberOfBars.  The  arc  connecting  the  two  events  is  a  scheduling  arc.  The  arc 
will  schedule  a  doBar(a)  event  some  time  units  after  the  execution  time  of  doFoo.  The 
time  of  scheduling  will  be  determined  by  the  distribution  of  tbar  written  next  to  the  arc 
with  no  parentheses.  If  no  time  distribution  is  written  next  to  the  arc,  the  scheduled  event 
will  be  scheduled  at  the  same  time  as  the  scheduling  event.  A  doBar(a)  event  will  only  be 
scheduled  if  the  condition  numberOfFoos  is  greater  than  numberOfBars,  is  met.  The  arc 
conditioning  is  written  in  parentheses  next  to  the  scheduling  arc.  the  event  doFoo  passes 
the  variable  “a”  to  the  event  “doBar(a),  written  in  rectangular  frame  next  to  the 
scheduling  arc. 

For  further  reading  on  SimKit  and  DES  see  [10]. 

2.  The  Simulation  Model 

The  following  two  figures  present  our  DES  event  graph.  Because  of  it  size,  the 
event  graph  is  divide  into  two  figures,  the  events  marked  in  red  are  common  to  both 
figures. 
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Figure  23 .  Simulation  event  graph  part  A 
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Figure  24.  Simulation  event  graph  part  B 
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C.  ANALYSIS 

1.  Positive  Salvo  Declaration  as  a  Function  of  the  Firing  Rule  ( M,T ) 

Our  MOE  is  the  proportion  of  positive  salvo  declarations  among  the  total  number 
of  salvos.  This  MOE  estimates  the  sensitivity  of  our  firing  rule  -  the  conditional 
probability  of  detecting  a  salvo  on  time.  We  call  this  MOE  the  probability  of  positive 
salvo  declaration. 

The  scenario  parameters  are: 

The  missile  is  “Dumb”.  There  are  N  =  5  rockets  in  the  salvo.  The  inter- firing  time 
is  uniformly  distributed  between  4.5  to  5.5  seconds.  The  time  to  escape  is  uniformly 
distributed  between  14  to  16  seconds;  The  probability  of  rocket  detection  is  a-  0.7;  to 
find  the  conditional  probability  of  kill  we  use  a  cookie  cutter  damage  function  with 
parameters  p  =  R2  =  1 .  We  assume  the  location  estimate  errors  are  iid  circular  normal 

r.v.  with  variance  cr2  =  1  .  The  time  between  salvos,  and  the  time  between  false  detections 
are  exponentially  distributed  with  mean  7200  seconds  (2  hours),  and  900  seconds  (15 
min)  respectively.  The  missile  flight  time  rm  =  20  seconds  is  fixed. 

Figure  25presents  the  probability  of  positive  salvo  declaration  for  different  fireing 

rules. 
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Probability  of  positive  salvo  decleration.  No  false  detections 


Figure  25.  Simulation  estimates*  of  the  Probability  of  positive  salvo  declaration  vs.  T 

(in  seconds)  for  different  values  of  M  ** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color 

First,  we  notice  that  the  probability  of  positive  salvo  declaration  is  increasing  with 
f  and  decreasing  with  M  .  This  is  the  expected  result;  as  we  increase  M  and  reduce  f 
we  make  our  firing  rule  tighter.  Take  for  instance  an  extreme  case  such  as(M  =  5,T  □  0) . 
This  is  the  tightest  firing  rule  in  Figure  25,  and  we  see  the  probability  of  positive  salvo 
declaration  is  zero.  On  the  other  hand,  (M  =  2,  T  =  40)  is  almost  the  least  tight  firing  rule 
in  the  figure.  We  see  the  probability  of  positive  salvo  declaration  is  almost  one. 
When(M  =  1) ,  the  probability  of  positive  salvo  declaration  is  of  course  one. 
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Finally,  we  observe  that  for  a  given  value  of  M  increasing  T  indefinitely 
(supposedly  making  the  firing  rule  infinitely  lenient),  does  not  lead  to  positive  salvo 
declaration  probability  of  one.  Instead  the  probability  increases  with  T  and  then  stays 
fixed  on  a  certain  value.  This  phenomenon  is  due  to  the  fact  that  in  counting  detections 
for  the  next  salvo  declaration,  we  discard  detections  used  for  previous  salvo  declarations. 
This  is  equivalent  to  saying  that  detections  that  generate  one  salvo  declaration  may  not 
overlap  with  detections  generating  the  next  salvo  declaration.  As  we  increase  T  we  do 
initially  make  the  firing  rule  less  demanding;  however  at  a  certain  point  it  becomes  more 
and  more  likely  that  the  time  window  will  include  detections  from  a  set  that  resulted  in  a 
previous  salvo  declaration.  From  this  point  on  increasing  T  will  not  make  any  difference 
since  any  detection  further  back  in  time  is  ignored  when  counting  detections  for  the  next 
salvo  declaration.  The  rate  of  positive  salvo  declarations  will  then  become  fixed 
regardless  of  any  increase  of  T  .  This  in  turn  leads  to  a  fixed  probability  of  positive  salvo 
declarations. 

2.  False  Salvo  Declaration  as  a  Function  of  the  Firing  Rule  ( M ,T ) 

The  MOE  for  the  false  declarations  is  the  ratio  between  the  number  of  false  salvo 
declarations  and  the  number  of  total  salvo  declarations  (including  duplicate  positive 
declarations  of  the  same  salvo).  We  call  this  MOE  the  false  declaration  proportion 
(FDP).  Another  MOE  is  the  false  declarations  rate  (FDR)  -  the  number  of  false 
declarations  per  unit  time  (FDR).  Figure  26and  Figure  27present  the  simulation  estimates 
for  the  FDP  and  the  overall  rate  of  false  salvo  declarations. 

The  scenario  parameters  are  as  in  0 
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FDP  for  diffrent  firing  rules 


Figure  26.  Simulation  estimates*  of  FDP  vs.  ( M  ,T )  ** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color 

The  area  circled  in  red  on  Figure  26shows  non-monotonic  behavior  and  missing 
points.  This  area,  having  the  tightest  firing  rules,  has  a  very  low  rate  of  salvo  declarations 
(false  or  positive).  The  low  rate  of  declarations  (see  next  figure  for  the  rate  of  false 
declarations)  makes  the  statistics  highly  variable,  and  in  some  cases  the  ratio  cannot  be 
computed  at  all  as  there  are  no  declarations  during  the  run.  This  makes  exploring  the 
region  of  low  false  declaration  probability  hard  to  do. 

Figure  27presents  an  MOE  which  may  be  more  useful  operationally  -  the  number 
of  false  salvo  declarations  per  unit  time  (Hour)  -  FDR.  This  translates  to  the  rate  of 
missiles  launched  unprovoked  per  hour. 
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4.5 


FDR  for  diffrent  firing  rules 


Figure  27.  Simulation  estimates*  FDR  vs.  ( M,T )  ** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color 

3.  Results  of  the  Base  Case  Scenario 

Dumb  missile 

Figure  28,  Figure  28and  Figure  29present  the  probability  of  success  the  FDR  and 
FDP  for  different  firing  rule  pairs  ( M ,T )  .  The  missile  is  “Dumb”  and  the  scenario  is  as 
in  Figure  25 
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The  probability  of  success  is  defined  here  as  in  Chapter  III.  The  simulation 
estimate  for  the  probability  of  success  is  the  ratio  between  the  number  of  salvos  the 
attacker  has  been  killed  and  the  total  number  of  salvos. 


Probability  of  success 


Figure  28.  Simulation  estimates*  for  the  Probability  of  success  vs.  T  (in  seconds)  for 

different  values  of  M  (“Dumb”  missile)** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color. 
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FDR  for  diffrent  firing  rules 


Figure  29.  Simulation  estimates*  of  FDR  vs.  T  (in  seconds)  for  different  values  of  M  ** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color. 
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FDP  for  diffrent  firing  rules 


Figure  30.  Simulation  estimates*  of  FDP  vs.  f  (in  seconds)  for  different  values  of  M  ** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color. 

For  clarity  of  presentation  only  the  cases  M  >2  are  presented  in  Figure  28.  We 
can  see  in  Figure  27  that  for  the  case  M  =  1  the  value  of  FDR  is  simply  the  rate  of  false 
detections  as  expected.  We  have  not  defined  the  defender’s  acceptable  value  of  FDR,  but 
it  seems  obvious  that  choosing  M  - 1  (effectively  responding  to  the  first  detection)  is 
unacceptable  due  to  the  high  value  of  FDR  (4 [FD/Hr])  ■  The  highest  probability  of  success 

(-0.55)  may  be  achieved  when  M  -2  however  when  M  =  2  there  is  still  a  significant 
FDR.  Since  the  FDR  is  monotone  increasing  with  T  ?  we  may  want  to  choose  T  slightly 
above  5  seconds.  This  will  give  the  lowest  value  of  FDR  possible  (0.02[/?D/&])  for  a  non 
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zero  probability  of  success  (-0.45).  Alternatively  the  defender  may  decide  to  settle  for 
lower  probability  of  success  and  choose  (M  =  3,  T  >  10) .  The  FDR  is  then  close  to  zero 
(<  0.005[ra/i?r])  and  the  probability  of  success  between  0.3  and  0.4. 

Method  1 :  Ground  Sensor  Guided  Missile 

The  following  figure  presents  the  simulation  estimates  of  the  probability  of 
success  for  the  ground  sensor  guided  missile.  The  rest  of  the  scenario  is  as  in  Figure  28 


Probability  of  success 


Figure  3 1 .  Simulation  estimates*  of  the  Probability  of  success  for  vs.  T  (in  seconds)  for 
different  values  of  M  (“Smart”  ground  sensor  guided  missile)** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color. 
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Changing  the  missile  type  does  not  influence  the  FDR  therefore  we  may  still  refer 
to  Figure  28and  Figure  27for  the  FDR.  The  most  significant  difference  from  Figure  28is 
in  the  curves  for  M  -  1  and  M  =  2  which  reach  considerably  higher  values.  However  this 
does  not  change  the  defender’s  decision  considerably;  the  defender  still  has  only  two 
viable  options  M  =  2  and  M  =  3 .  The  defender  may  be  more  inclined  to  choose  M  -  2 
because  of  the  increase  in  the  probability  of  success. 

Method  2:  Onboard  sensor  Guided  Missile 

The  following  figure  presents  the  probability  of  success  for  the  onboard  sensor 
guided  missile;  the  sensor  is  “Passive”  and  receives  10  times  less  energy  than  the  ground 
sensor.  The  rest  of  the  scenario  is  as  in  Figure  28 


Probability  of  success 


Figure  32.  Simulation  estimates*  of  Probability  of  success  for  vs.  T  for  different  values 

of  M  (“Smart”  onboard  “Passive”  sensor  guided  missile)** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 


rocket  salvos. 
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**To  properly  read  this  figure,  print  in  color. 

Changing  the  missile  type  does  not  influence  the  FDP  therefore  we  may  still  refer 
to  Figure  28and  Figure  27for  the  FDP.  Changing  from  the  “dumb”  missile  to  an  onboard 
“Passive”  sensor  missile  does  not  affect  the  defender’s  decision  at  all.  It  seems  the  only 
difference  between  Figure  28and  Figure  32is  in  the  curve  for  M  =  1  which  is  not  an 
option  because  of  the  high  FDR. 

It  seems  given  only  these  three  options  that  the  defender  would  choose  to  use  the 
ground  sensor  guided  missile,  possibly  using  the  firing  rule  pair  (M  =  3,  f  >  10) . 

4.  The  Effect  of  the  False  Detection  Rate 

No  False  detections 

Let  us  first  look  at  the  extreme  case  where  the  rate  of  false  detections  is  zero.  The 
following  figure  presents  the  simulation  estimates  for  the  probability  of  positive  salvo 
declaration  when  there  are  no  false  detections  (False  detections  are  simply  not  scheduled 
in  these  simulation  runs).  Although  there  are  no  false  detections  the  defender  still 
implements  the  same  firing  rule.  The  rest  of  the  scenario  is  the  same  as  in  Figure  28  Each 
estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000  rocket  salvos. 
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Probability  of  positive  satvo  decleration.  No  false  detections 


Figure  33.  Simulation  estimates*  of  Probability  of  positive  salvo  declaration  vs.  T  for 

different  values  of  M  No  false  detections  ** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color. 

When  M  -  1 ,  the  probability  of  positive  salvo  declaration  is  always  1.  For  M  >  1 
we  notice  there  are  sudden  increases  in  the  probability  as  f  is  increased.  Take  for 
example  the  case  where  M  -  2 .  Recall  the  inter-firing  time  is  uniformly  distributed 

between  4.5  and  5.5  seconds.  Therefore  when  the  decision  variable  T  <4.5  seconds, 
there  is  no  combination  of  detected  rockets  within  the  salvo  that  can  generate  a  salvo 
declaration.  Any  two  detected  rockets  will  be  detected  at  a  greater  time  interval  than  T . 
However  once  T  >  5.5 ,  every  two  consecutive  detected  rockets  will  fall  within  the  time 
interval  T  and  therefore  will  trigger  a  salvo  declaration.  Notice  that  the  first  increase  in 
the  probability  of  salvo  declaration  for  M  -  2  is  about  equal  to  the  probability  of 
detecting  any  combination  containing  at  least  two  consecutive  detected  rockets,  given 
by:  P  =  l-(l-a)5  -5  -a  -(1-a)5  -6 -a2  -(1-a)3  -a3  -(1-a)2  =  0.8859 .  As  we  increase 


64 


f  more  combinations  of  detected  rockets  may  generate  a  salvo  declaration  and  so  the 
probability  of  salvo  declaration  increases.  The  maximal  probability  of  salvo  declaration 

for  M -2  (given  when  T>  20)  is  about  .997  which  corresponds  to  one  minus  the 
probability  of  that  no  rocket  is  detected  in  the  salvo  given  by:  P  =  l-(l-a)5  =  0.9976  , 

the  only  combination  which  will  not  trigger  a  salvo  declaration  when  M  =  2  ,  T  >  20 . 
High  Rate  of  False  Detections 

Let  us  now  explore  the  effect  of  increasing  the  ground  sensor’s  rate  of  false 
detections.  The  following  figure  presents  the  probability  of  success  and  the  FDP.  The 
mean  time  between  false  detection  is  60  seconds.  The  rest  of  the  scenario  is  the  same  as 
in  Figure  28 


Probability  of  success 
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Figure  34.  Simulation  estimates*  for  the  Probability  of  success  vs.  T  (in  seconds)  for 
different  values  of  M  (“Dumb”  missile).High  rate  of  false  detections  ** 


65 


Probability  of  false  salvo  (federation 


Figure  35.  Simulation  estimates*  of  FDP  vs.  T  (in  seconds)  for  different  values  of  M  . 

High  rate  of  false  detections  ** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  these  figures,  print  in  color. 

Increasing  the  rate  of  false  detections  increases  the  FDP  (compare  to  Figure  26). 
False  detections  are  the  cause  for  false  salvo  declaration  and  so  this  result  is  expected. 
Comparing  Figure  34to  Figure  28  we  see  that  the  probability  of  success  is  increased  as 

well  (this  is  most  obvious  when  M  =  5).  To  understand  the  contribution  of  false 
detections  to  the  probability  of  success,  think  of  the  case  when  M  =  4  .  When  M  =  4  and 
in  the  absence  of  false  detections,  the  firing  rule  demands  for  4  rockets  to  be  detected. 
The  probability  of  this  event  is  not  very  high 
(N-a4  -{\-a)  +  a5  =  5-0.74  -0.3  +  0.75  =0.52).  Furthermore  in  the  uniform  distribution 
case,  if  the  decision  variable  T  is  low  enough  there  will  never  be  enough  time  for  all  4 
rockets  to  be  fired  within  T .  This  reduces  further  the  probability  of  salvo  declaration. 
Indeed  we  see  that  in  figure  Figure  28whenM  =  4  and  f  <  (M  -1)  ■  ju  .  =  1 5  then  the 
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probability  of  success  is  close  to  zero.  However  in  the  presence  of  false  detections,  fewer 
than  4  rockets  in  the  salvo  need  to  be  detected  for  a  salvo  to  be  declared.  Therefore 
increasing  the  rate  of  false  detections  increases  the  probability  of  salvo  declaration  and  in 

turn  the  probability  of  success.  Notice  for  M  =  4 and  10<r<15  in  Figure  34  the 
probability  of  success  is  significant  and  not  close  to  zero  as  in  Figure  28  We  learn  from 
this  discussion  that  the  probability  of  success  presented  in  Chapter  III  with  no  false 
detections,  gives  a  lower  bound  on  the  maximal  probability  of  success  over  T  for  the 
case  with  false  detections. 

FDP  has  sudden  drops  at  some  values  of  T .  These  drops  correspond  to  sudden 
increases  in  the  probability  of  success.  As  we  have  seen  earlier  increasing  T  allows  more 
combinations  of  detected  rockets  in  the  salvo  to  generate  a  salvo  declaration.  This 
suddenly  increases  the  rate  of  positive  salvo  declaration  but  only  slightly  changes  the  rate 
of  false  salvo  declarations.  The  net  result  is  a  sudden  drop  in  the  FDP. 

The  “wild”  swings  in  the  FDP  for  low  values  of  T  when  M  -  5  are  due  to 
simulation  variability  as  explained  in  4.C.2. 

5.  The  Effect  of  the  Missile  Flight  Time 

The  response  of  the  probability  of  success  to  variations  in  the  missile  flight  time 
has  been  shown  in  Chapter  III  in  the  case  with  no  false  detections.  We  have  already 
concluded  that  the  results  in  Chapter  III  give  a  lower  bound  on  the  maximal  probability 
of  success  (over  T )  with  false  detections.  For  relatively  low  rate  of  false  detections,  we 
expect  the  maximal  values  of  the  probability  of  success  to  be  close  to  the  ones  in  the  case 
with  no  false  detections.  For  instance,  in  Figure  28the  maximal  probability  of  success 

when  M  =  2  is  about  0.55.  Comparison  of  this  value  to  the  probability  of  success  when 
M  =  2  and  rm  =  20  in  Figure  7shows  they  are  about  the  same.  We  can  repeat  this  exercise 

for  any  value  of  M  and  for  any  type  of  missile  (comparing  Figure  llto  Figure  31  and, 
Figure  19  to  Figure  32)  We  demonstrate  this  one  more  time  by  changing  the  missile  flight 
time  to  rm  =  2  and  again  comparing  with  Figure  7 
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Probability  of  success 


Figure  36.  Simulation  estimates*  of  Probability  of  success  for  vs.  T  for  different  values 

of  M  (“Dumb”  missile).  rm  =  2  ** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  this  figure,  print  in  color. 

6.  The  Effect  of  the  Salvos  Rate 

Let  us  now  explore  the  effect  of  increasing  the  rate  of  salvos.  The  following  two 
figures  present  estimates  of  the  probability  of  success  and  the  FDP.  Each  estimate  is  from 
a  single  simulation  run  containing  a  sequence  of  1000  rocket  salvos.  The  mean  time 
between  salvos  is  900  seconds  (15  minutes),  a  substantial  decrease  from  2  hours.  The  rest 
of  the  scenario  is  the  same  as  in  Figure  28 

Figure  37and  Figure  28are  the  same.  Increasing  the  salvo  rate  does  not  change  the 
probability  of  success.  Comparing  Figure  37and  Figure  29we  notice  the  FDP  drops 
dramatically.  Recall  the  definition  of  a  false  salvo  declaration  basically  says  a  false  salvo 
declaration  is  a  salvo  declaration  which  occurs  during  the  time  between  salvos. 
Decreasing  the  mean  time  between  salvos  means  that  for  a  given  time  interval  during  the 
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simulation  run,  there  would  be  on  average  more  salvos  and  less  time  during  which  false 
salvos  may  occur.  Since  the  firing  rule  is  unchanged  the  rate  of  positive  salvo 
declarations  during  the  salvo  and  the  rate  of  false  salvo  declaration  between  salvos  are 
unchanged.  Overall  the  simulation  run  will  have  more  positive  salvo  declarations  and  less 
false  salvo  declarations.  The  net  result  is  a  lower  FDP.  It  is  interesting  to  note  that 
although  the  FDR  has  not  changed  the  FDP  has  changed.  If  the  defender  is  more 
concerned  with  the  ratio  of  false  to  positive  salvo  declarations  rather  than  the  rate  of  false 
salvo  declarations  per  unit  time,  he  may  choose  a  less  conservative  firing  rule  under  a 
more  intensive  rocket  attack.  In  other  words  when  the  attacker  shoots  more  salvos,  the 
defender  may  become  more  “happy  trigger”.  This  last  observation  means  that  when  the 
rate  of  salvo  is  higher,  the  defender  may  choose  a  more  lenient  firing  rule.  Perhaps  in  the 
case  above  the  defender  will  choose  (M  =  2 ,f  =  6) . 
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Probability  of  success 


Figure  37.  Simulation  estimates*  of  the  Probability  of  success  for  vs.  T  (in  seconds)  for 

different  values  of  M  (“Dumb”  missile)  High  salvo  rate** 


Probability  of  false  salvo  decleration 


Figure  38.  Simulation  estimates*  of  the  FDP  vs.  T  (in  seconds)  for  different  values  of 

M  High  salvo  rate  ** 
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*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  these  figures,  print  in  color. 

7.  The  Exponential  Case 

We  now  explore  the  case  where  the  inter- firing  time  and  escape  time  are 
exponentially  distributed.  The  following  two  figures  present  the  probability  of  success 
and  the  FDP.  The  inter-firing  time  and  escape  time  are  exponentially  distributed  with 
means  5  and  15  seconds  respectively.  The  rest  of  the  scenario  is  the  same  as  in  Figure  28 

Comparing  the  exponential  case  to  the  uniform  case  (Figure  28and  Figure  28),  we 
notice  how  in  the  exponential  case  the  curves  become  smoother  and  have  no  sharp  drops 
or  increases  as  in  the  unifonn  case.  The  FDP  does  not  present  non  monotone  behavior 
when  rm  is  very  low.  In  the  exponential  case  the  probability  of  success  is  never  zero  no 
matter  how  demanding  the  firing  rule  is.  Comparing  Figure  39to  Figure  9we  can  see 
again  how  the  maximal  probability  of  success  over  T  corresponds  to  the  values  of  the 
probability  of  success  with  no  false  detections. 
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Probability  of  success 


Figure  39. 


Simulation  estimates*  of  the  Probability  of  success  for  vs.  T  for  different 
values  of  M  (“Dumb”  missile)  Exponential  case** 
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Probability  of  false  salvo  decleration 


Figure  40.  Simulation  estimates*  of  FDP  vs.  T  for  different  values  of  M 

Exponential  case  ** 

*  Each  estimate  is  from  a  single  simulation  run  containing  a  sequence  of  1000 
rocket  salvos. 

**To  properly  read  these  figures,  print  in  color. 
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V.  SUMMARY 


In  this  thesis,  we  have  explored  and  discussed  tactical  considerations  in 
countering  mortar  and  rocket  fire,  when  the  enemy  has  adopted  a  “shoot  and  scoot” 
tactic. 

The  probability  of  success  -  the  probability  of  a  timely  effective  hit  of  the  attacker 
-  is  used  as  the  MOE  for  the  defender’s  effectiveness.  We  have  used  the  probability  of 
false  salvo  declaration  proportion  (FDP)  and  the  false  salvo  declarations  rate  (FDR)  as 
MOEs  for  the  defender’s  chances  of  mistakenly  firing  a  missile  and  increasing  the 
collateral  damage. 

We  have  developed  a  tool  for  comparing  different  defender  tactics  and  responses 
(firing  rules).  In  Chapter  III,  we  have  used  an  analytic  model  under  the  assumption  that 
there  are  no  false  detections  (and  hence  zero  FDP  and  FDR).  This  analytic  model  allowed 
us  to  quantify  the  defenders  dilemma  between  waiting,  gathering  more  information  and 
increasing  accuracy,  and  responding  early  at  the  expense  of  more  information  and  with 
less  accuracy. 

The  analytic  model  presented  how  the  probability  of  success  depends  on  the 
sensor’s  probability  of  rocket  detection,  the  missile  flight  time  and  the  attacker’s  inter¬ 
firing  and  escape  time  distributions. 

We  have  explored  a  suggested  new  tactic  where  the  defender  launches  his  missile 
as  soon  as  possible  based  on  a  rough  estimate  of  the  attacker’s  location.  The  missile’s 
trajectory  is  later  corrected  and  updated  while  the  missile  is  airborne,  using  ground  sensor 
guidance  or  different  types  of  onboard  sensors.  The  effectiveness  of  the  suggested  tactic 
depends  on  the  parameters  of  the  scenario;  however  it  is  most  useful  when  the  missile 
flight  time  is  long  compared  to  the  mean  salvo  length  and  escape  time. 

In  Chapter  IV,  we  have  introduced  false  detections  to  the  model.  To  deal  with  this 
predicament  we  have  built  a  simulation.  The  defender’s  firing  rule  has  been  extended  to 
account  for  false  detections  by  introducing  a  sliding  time  window.  The  problem  in 
Chapter  IV  was  to  find  a  firing  rule  which  gives  the  highest  possible  probability  of 
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success  while  keeping  FDP  and  FDR  as  low  as  possible.  It  surprising  to  note  that 
presence  of  false  detections,  actually  contributes  the  defender’s  probability  of  success. 
This  happens  because  false  detections  trigger  the  defender  to  launch  his  missile  in  cases 
he  otherwise  wouldn’t  have,  even  though  a  salvo  is  ongoing  (e.g.  no  rockets  were 
detected).  Of  course  this  improvement  comes  at  the  cost  of  mistakenly  tiring  missiles 
causing  collateral  damage. 

We  discussed  the  simulation  results  for  a  base  case  and  varied  the  rate  of  salvos 
the  rate  of  false  detections  and  the  missile  flight  time  for  sensitivity  analysis. 

Follow-on  research  may  include  taking  the  salvo  rate,  inter- firing  and  escape  time 
distributions  from  real  life  data  and  exploring  the  ideas  presented  in  this  paper  on  that 
data.  It  would  also  be  interesting  to  survey  existing  weapons  and  sensors  and  their 
parameters  to  see  if  it  possible  to  use  any  of  them  for  the  new  tactic  suggested  in  this 
thesis. 
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APPENDIX  A 


Although  the  results  for  the  Normal  case  are  very  similar  to  the  results  of  the 
Uniform  case,  we  present  here  for  the  purpose  of  completeness  our  results  for  the  case 
where  the  inter-firing  and  escape  time  are  Normally  distributed. 

There  are  N  =  5  rockets  in  the  salvo.  The  inter- firing  time  is  normally  distributed 

1 

with  mean  5  seconds  and  variance  —sec".  The  time  to  escape  is  normally  distributed 

1 

with  mean  15  seconds  and  variance  -see  .  Notice  that  the  inter  firing  and  escape  time 
variances  and  means  are  the  same  as  in  Figure  IThe  rest  of  the  scenario  is  as  in  Figure  1 


Probability  of  success  for  different  missile  flight  times 


Figure  4 1 .  Probability  of  success  for  different  missile  flight  times 
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Probability  of  timely  hit  for  different  missile  flight  times  and  different  number  of  rockets  left 


Figure  42.  Probability  of  a  timely  hit  for  different  missile  flight  times  and  number  of 

residual  rockets* 

*  To  properly  read  this  figure,  print  in  color. 

As  mentioned  above  comparing  Figure  land  Figure  2with  Figure  4 land  Figure 
42We  see  that  the  results  for  the  Normal  case  are  very  similar  to  the  results  of  the 
Uniform  case. 
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